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On Nonlinear Asymptotic Stability of the Lane-Emden Solutions for the 

Viscous Gaseous Star Problem 

Tao Luo, Zhouping Xin, Huihui Zeng 

Abstract 

This paper proves the nonlinear asymptotic stability of the Lane-Emden solutions for spher¬ 
ically symmetric motions of viscous gaseous stars if the adiabatic constant 7 lies in the stability 
range (4/3,2). It is shown that for small perturbations of a Lane-Emden solution with same 
mass, there exists a unique global (in time) strong solution to the vacuum free boundary prob¬ 
lem of the compressible Navier-Stokes-Poisson system with spherical symmetry for viscous stars, 
and the solution captures the precise physical behavior that the sound speed is C 1 ' 2 -Holder con¬ 
tinuous across the vacuum boundary provided that 7 lies in (4/3, 2). The key is to establish the 
global-in-time regularity uniformly up to the vacuum boundary, which ensures the large time 
asymptotic uniform convergence of the evolving vacuum boundary, density and velocity to those 
of the Lane-Emden solution with detailed convergence rates, and detailed large time behaviors 
of solutions near the vacuum boundary. In particular, it is shown that every spherical surface 
moving with the fluid converges to the sphere enclosing the same mass inside the domain of 
the Lane-Emden solution with a uniform convergence rate and the large time asymptotic states 
for the vacuum free boundary problem ( 11 . 1 . 21 ) are determined by the initial mass distribution 
and the total mass. To overcome the difficulty caused by the degeneracy and singular behavior 
near the vacuum free boundary and coordinates singularity at the symmetry center, the main 
ingredients of the analysis consist of combinations of some new weighted nonlinear function¬ 
als (involving both lower-order and higher-order derivatives) and space-time weighted energy 
estimates. The constructions of these weighted nonlinear functionals and space-time weights 
depend crucially on the structures of the Lane-Emden solution, the balance of pressure and 
gravitation, and the dissipation. Finally, the uniform boundedness of the acceleration of the 
vacuum boundary is also proved. 
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1 Introduction 


1.1 Problem 


In the fundamental hydrodynamical setting (cf. [I]), the evolving boundary of a viscous gaseous 
star (the interface of fluids and vacuum states) can be modeled by the following free boundary 
problem of the compressible Navier-Stokes-Poisson equations: 


p t + div(pu) = 0 

(pu) i + div(pu <g> u) + div© 

p > 0 

p = 0 and ©n = 0 
V(r (i)) = u-n, 

(p,u) = Oo,u 0 ) 


-pVx'k 


in 

in £l(t), 
in £l(t), 

on T (t) := d£l(t), 


on O := 0 ( 0 ). 


( 1 - 1 - 1 ) 


Here (x, t) € M 3 x [0, oo), p, u, & and 'F denote, respectively, the space and time variable, density, 
velocity, stress tensor and gravitational potential; O(t) C M 3 , r(f), V(r(f)) and n represent, respec¬ 
tively, the changing volume occupied by a fluid at time t, moving interface of fluids and vacuum 
states, normal velocity of T(t) and exterior unit normal vector to T(t). The gravitational potential 
is described by 

\F(x, t) = — G [ r-^—^dy satisfying = AirGp in O(t) 

Jn(t) l x - y| 

with the gravitational constant G taken to be unity for convenience. The stress tensor is given by 


6 = ph 


Ai 


^Vu + Vu* 


2 

3 


(divu )I 3 


A 2 (divu )/ 3 , 


where 1$ is the 3x3 identical matrix, p is the pressure of the gas, Ai > 0 is the shear viscosity, 
A 2 > 0 is the bulk viscosity, and Vu f denotes the transpose of Vu. We consider the polytropic, 
gases for which the equation of state is given by 


P = P(P) = KP 1 , 

where K > 0 is a constant set to be unity for convenience, 7 > 1 is the adiabatic exponent. 

For a non-rotating gaseous star, it is important to consider spherically symmetric motions since 
the stable equilibrium configurations, which minimize the energy among all possible configurations 
(cf. [2lj), are spherically symmetric, called Lane-Emden solutions. In this work, we are concerned 
with the three-dimensional spherically symmetric solutions to the free boundary problem (jl.l.ip and 
its nonlinear asymptotic stability toward the Lane-Emden solutions. The aim is to prove the global- 
in-time regularity uniformly up to the vacuum boundary of solutions when 4/3 < 7 < 2 (the stable 
index) capturing an interesting behavior called the physical vacuum (cf. [3l 0J [GO, ,15; 231 [Ml 421) 
which states that the sound speed c = \Jp'(p) is C 1//2 -Holder continuous near the vacuum boundary, 
as long as the initial datum is a suitably small perturbation of the Lane-Emden solution with the 
same total mass. Furthermore, we establish the large time asymptotic convergence of the global 
strong solution, in particular, the convergence of the vacuum boundary and the density, to the the 
Lane-Emden solutions with the detailed convergence rate as the time goes to infinity. 
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In the spherically symmetric setting, that is, £l(t) is a ball with the changing radius R(t), 
p(x, t) = p(r, t) and u(x, t) = u(r, t)~x./r with r = |x| G (0, R(t )); 
system ( 11 . 1 . 11 ) can then be rewritten as 


to 

+ 

to 

■i 

II 

O 

in 

(0, R(t)), 

( 1 . 1 . 2 a) 

f r (( r2u )r\ 

p(u t + uu r ) + Pr + 47T pr 2 J p(s , t)s' 2 ds = p, ( J 

in 

(0, R(t)), 

( 1 . 1 . 2 b) 

p > 0 

in 

[0, R(t)), 

( 1 . 1 . 2 c) 

4 / u \ / u\ 

p = 0 and -Ai yu r -J + A 2 yu r + 2 —J = 0 

for 

ho' 

II 

( 1 . 1 . 2 d) 

R(t) = u(R(t),t ) with R( 0) = Ro, u(0,t) = 0, 



( 1 . 1 . 2 e) 

(p,u) = {po,u 0 ) 

on 

(0, Ro), 

( 1 . 1 . 2 f) 


where p = 4Ai/3 + A 2 > 0 is the viscosity constant. (I1.1.2cl) and ( 11 . 1 . 2 dl) state that r = R(t) is the 
vacuum free boundary at which the normal stress 6 n = 0 reduces to 

4 / u \ / u \ 

V — gAi ya r -J — A 2 yu r + 2—j =0 for r = R(t), t > 0; 

(| 1 . 1 . 2 cl) describes that the free boundary issues from r = Ro and moves with the fluid velocity, and 
the center of the symmetry does not move. The initial domain is taken to be a ball {0 < r < Ro}, 
and the initial density is assumed to satisfy the following condition: 

po{r) >0 for 0 < r < Ro, po{Ro ) = 0 and — 00 < < 0 at r = Ro ; (1.1.3) 

so 


Po 1 ( r ) ~ Ro — r as r close to Ro, (1.1.4) 

that is, the initial sound speed is C 2 -Holder continuous across the vacuum boundary. The unknowns 
here are p, u and R(t). 

The requirement (11.1.31) for the initial density near the vacuum boundary is motivated by that 
of the Lane-Emden solution, p, (cf. [l ;! [22j) which solves 

pr 

d r (fP) + ^Trr~ 2 p / p(s)s 2 ds = 0. (1.1.5) 

Jo 

The solutions to (jl.1.51) can be characterized by the values of 7 (cf. [[22]) for given finite total mass 
M > 0, if 7 £ (6/5,2), there exists at least one compactly supported solution. For 7 G (4/3,2), 
every solution is compactly supported and unique. If 7 = 6/5, the unique solution admits an 
explicit expression, and it has infinite support. On the other hand, for 7 G (1,6/5), there are no 
solutions with finite total mass. For 7 > 6/5, let R be the radius of the stationary star giving by 
the Lane-Emden solution, then it holds (cf. [|22( 131]) 

p 7_1 (r) ~ R — r as r close to R. (1.1.6) 
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1.2 Motivations and goals 

The problem of nonlinear asymptotic stability of Lane-Emden solutions is of fundamental impor¬ 
tance in both astrophysics and the theory of nonlinear PDEs. It is believed by astrophysicists 
that Lane-Emden solutions are stable for 4/3 < 7 < 2 (cf. pD, 39]) since they minimize the total 
energy among all the possible configurations. The main aim of this paper is to justify rigorously 
the precise sense of this stability. In fact, we prove for the viscous gaseous star with 4/3 < 7 < 2 , 
the Lane-Emden solution is strongly stable in the sense that it is asymptotically nonlinear stable. 
The first step for this purpose is to prove the global existence of strong solutions. However, due to 
the high degeneracy of system (11.1.21) caused by the behavior (ll.l.3|) near the vacuum boundary, it 
is a very challenging problem even for the local-in-time existence theory. Indeed, the local-in-time 
well-posedness of smooth solutions to vacuum free boundary problems with the behavior that the 
sound speed is C' 1 / 2 -Holder continuous across vacuum boundaries was only established recently 
for compressible inviscid flows (cf. [3] @1 ITU ll~5j ) (see also [29] for a local-in-time well-posedness 
theory in a new functional space for the three-dimensional compressible Euler-Poisson equations in 
spherically symmetric motions). For the vacuum free boundary problem ([1.1.21) of the compress¬ 
ible Navier-Stokes-Poisson equations featuring the behavior (|1.1.3I) near the vacuum boundary, a 
local-in-time well-posedness theory of strong solutions was established in [12] . In order to obtain 
the nonlinear asymptotic stability of Lane-Emden solutions, it turns out that suitable estimates for 
higher order derivatives uniformly up to the vacuum boundary are necessary. Indeed, this turns 
out to be essential to prove the convergence of the evolving vacuum boundary and the uniform con¬ 
vergence of the density to those of Lane-Emden solutions, in addition to the uniform convergence 
of the velocity. We show the global-in-time regularity of solutions when 4/3 < 7 < 2 capturing 
the behavior (11.1.41) (or ([1,1.61) ) when the initial data are small perturbations of and have the same 
total mass as the stationary solution, p, given by (11.1.51) . It should be remarked that the regularity 
estimates near boundaries are notoriously difficult. This is particularly so for the vacuum boundary 
problem (II.1.21) due to the high degeneracy caused by the singular behavior of ([1.1.3D near vacuum 
states. 

Our nonlinear asymptotic stability results can be stated more precisely as follows. Suppose 
that the initial datum (po,Uo, Rq) is a small perturbation of the Lane-Emden solution (p,0,R) in 
a suitable sense (see Theorem 12. 3 1) and has the same total mass, 

rRo rR 

/ r 2 po(r)dr = / r 2 p(r)dr, 

Jo Jo 

then there is a unique global-in-time strong solution (p, u , R(t)) (0 < t < + 00 ) to ([1.1.2D which is 
regular uniformly up to the vacuum boundary r = R(t). Moreover, let r(x,t ) be the radius of the 
ball inside B R ^ (0) satisfying: 

r t (x,t) = u(r(x,t),t) and r(x, 0) = ro(a;) for 0 < x < R, (1.2.1) 



Then 



for 0 < x < R. 


lim 

t—yoo 


{r{x,t)-x, p(r(x,t),t) -p(x), u(r(x, t), t)) || l - ([0 ,r] ) = 0 


( 1 . 2 . 2 ) 


(1.2.3) 


with some detailed convergence rates. Notice that (11.2.11) means that the sphere r = r(x, t ) with the 
initial position r = r q(x) is moving with the fluid and (11.2.21) means that the initial mass inside the 
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ball B ro ( x }( 0) is the same as that of the Lane-Emden solution inside the ball B x { 0) for 0 < x < R. 
It follows from the conservation of mass that the mass inside the ball B r i xt \{ 0) at the instant t is 
the same as that of the Lane-Emden solution inside the ball B x { 0) for 0 < x < R. In particular, 
the vacuum boundary is given by 

R(t) = r(R, t). 

The convergence of r(x,t ) to x in ([1.2.31) means that every spherical surface moving with the fluid 
converges to that inside the domain of Lane-Emden solution enclosing the same mass, in particular, 
the evolving vacuum boundary R(t) converges to the vacuum boundary R as time goes to infinity. 
This also gives the large time asymptotic convergence of every particle moving with the fluid since 
the motion is radial. Moreover, the convergence (11.2.31) means that the large time asymptotic states 
for the free boundary problem (11.1.2(1 are determined completely by the initial mass distribution 
and total mass. Besides the convergence (11.2.31) . we also establish convergence rates of higher- 
order norms involving derivatives, and show that the vacuum boundary R(t) has the regularity 
of W 2 ’°°([0, oo)) under a compatibility condition of the initial data with the boundary condition 
which implies that the acceleration of the vacuum boundary is uniformly bounded for t G [0, oo). 
(Indeed, one may check from the proof that every particle moving with the fluid has the bounded 
acceleration for t G [0, oo).) These results give a rather clear and complete characterization of the 
behavior of solutions both in large time and near the vacuum boundary. 

The results obtained in the present work are among few results of global strong solutions to 
vacuum free boundary problems of compressible fluids capturing the singular behavior of (11.1. 3D . 
which is difficult and challenging due to the degeneracy caused by the physical vacuum and co¬ 
ordinates singularity at the center of the symmetry. We overcome this difficulty by establishing 
higher-order estimates involving the second-order derivatives of the velocity field, together with 
decay estimates of lower-order norms. This is achieved by combining some new weighted nonlinear 
functionals (involving both lower-order and higher-order derivatives) and space-time weighted en¬ 
ergy estimates. The constructions of these weighted nonlinear functionals and space-time weights 
depend crucially on the structure of Lane-Emden solutions (in particular, the behavior (11.1.61) near 
the vacuum boundary), the balance between the pressure and self-gravitation, and the dissipation 
of the viscosity. 

We sketch here the main ideas and methods used in this article. The original free boundary 
problem (11.1.21) is reduced to an initial boundary value problem on a fixed domain x G [0, .R] by the 
Lagrangian particle trajectory formulation (II. 2. ID and (II. 2. 2D with R being the radius of the Lane- 
Emden solution so that the domain of the Lane-Emden solution becomes the reference domain. In 
this formulation, essentially the basic unknown is the particle trajectory r(x,t ) defined by (II. 2. ID 
and (II. 2. 2D (more precisely, the radius of each evolving surface inside the evolving domain, which 
is called the particle trajectory for simplicity here and from now on), by which the density and 
velocity are determined. For problem (jl.l.2p . this formulation is preferred because one can use 
it to trace each particle in the evolving domain, in particular, the evolving vacuum boundary. 
(Indeed, the approach of using the Lagrangian particle trajectory formulation, i.e., using the flow 
map of the Eulerian velocity field, to reduce free boundary problems of compressible fluids set on 
time-dependent domains to fixed time-independent domains was first adopted in [3j, and latter on 
in urns in]-) To make this strategy work, a crucial point is to obtain the positive lower and 
upper bounds of the derivative of the particle trajectory, which are derived by a pointwise estimate 
away from the center of symmetry, and an interior L 2 -estimate of its second derivative. Various 
multipliers are applied to establish the decay estimates of lower-order norms and the regularity near 
the vacuum boundary. In order to obtain the higher-order estimates, we first study the problem 
obtained by differentiating the original problem in the tangential direction, and then get weighted 
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estimates for the second derivatives of the velocity field using the viscosity term. However, due to 
the degeneracy of (11.1.311 , the dissipation of the viscosity alone is not enough for the global-in-time 
estimates, and we have to make full use of the balance between the pressure and gravitation. For 
this purpose, we decompose the gradient of the pressure as two parts, the first part is to balance the 
gravitation and the second part is an anti-derivative of the viscosity along the particle trajectory 
with a degenerate weight to match the viscosity. The main ideas and strategy of establishing the 
decay estimates and high-order regularity estimates will be given in Section 13.21 

1.3 Review of related works 

There have been extensive works on the studies of the Euler-Poisson and the Navier-Stoke-Poisson 
equations with vacuum, especially in recent years. We will concentrate on those closely related to 
the stability of vacuum dynamics. The stability problem has been important in the theory of gaseous 
stars which has been studied extensively by astrophysicists (cf. [1, Hi (T9]). The linear stability 
of Lane-Emden solutions was studied in [22]. A conditional nonlinear Lyapunov type stability 
theory of stationary solutions for 7 > 4/3 was established in [36] using a variational approach, 
by assuming the existence of global solutions of the Cauchy problem for the three-dimensional 
compressible Euler-Poisson equations (the same type of nonlinear stability results for rotating stars 
were given by EH m- For 7 £ (6/5,4/3), the nonlinear dynamical instability of Lane-Emden 
solutions was proved by m and m in the framework of free boundary problems for Euler-Poisson 
systems and Navier-Stokes-Poisson equations, respectively. A nonlinear instability for 7 = 6/5 was 
proved by m ■ For 7 = 4/3, an instability was identified in [5] that a small perturbation can cause 
part of the mass to go off to infinity for inviscid flows. 

It should be noted that the stability result in [36] is in the framework of initial value problems in 
the entire M 3 -space and involves only a Lyapunov functional which is essentially equivalent to a L p - 
norrn of difference of solutions, and the vacuum boundary cannot be traced. Another interesting 
work is on the vacuum free boundary problem of modified compressible Navier-Stokes-Poisson 
equations with spherical symmetry (cf. [TO]), where the existence of a global weak solution was 
proved for a reduced initial boundary value problem after using the Lagrangian mass coordinates, 
under some constraints on the ratio of the coefficients of the shear viscosity and bulk viscosity. 
In contrast to the strong stabilitv result in (11.2.31). for the global weak solutions obtained in m, 
only the uniform convergence of the velocity u(r,t ) is proved, due to the lack of regularity near 
the vacuum boundary. The ideas and techniques developed in this paper can be applied to this 
modified compressible Navier-Stokes-Poisson equations to obtain a strong stability result as in 
([1.2.31) . Indeed, our global-in-time regularity gives not only the decay estimates for the weighted 
norms \\p' y ^ 2 (r x (x,t) — l)||/y ([ 0) fi]) and H^'CxII/y ([o,r]) as in [TO] , but also the decay estimates of the 
unweighted norms of || r x (x,t) — 1||z,2([o,.r])> H u x||l°°([o,.r]) an d some uniform estimates on the second 
derivatives valid up to the vacuum boundary, which are crucial to the nonlinear asymptotic stability 
for this modified model. Furthermore, our theory holds without the restrictions on the viscosity 
coefficients as in m- This will be reported in a forthcoming paper (cf. [30]). 

The results obtained in the present work are for spherically symmetric motions. For the problem 
of general three dimensional perturbations, in addition to motions in the normal direction, one will 
have to estimate motions in the tangential direction. Besides this, the evolution of the geometry 
of the free surface will be estimated also. We believe the ideas and techniques developed in this 
paper will be useful to the study of the problem of general three dimensional perturbations. 

We conclude the introduction by noting that there are also other prior results on free boundary 
problems involving vacuum for compressible Navier-Stokes equations besides the ones aforemen¬ 
tioned. For the one-dimensional motions, there are many results concerning global weak solutions 
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to free boundary problems of the Navier-Stokes equations, one may refer to [si E5i Eai da e da 
mmm and references therein. As for the spherically symmetric motions, global existence and 
stability of weak solutions were obtained in [33], El] to compressible Navier-Stokes equations for 
gases surrounding a solid ball (a hard core) without self-gravitation (see also [ 2 ] for a compressible 
heat-conducting flow). However, those results are restricted to cut-off domains excluding a neigh¬ 
borhood of the origin. It should be noted that for a modified system of Navier-Stokes equations, a 
global existence of weak solutions with spherical symmetry containing the origin was established in 
du for which the density does not vanish on the boundary. For a class of free boundary problems 
of compressible Navier-Stokes-Poisson equations away from vacuum states, the readers may refer 
to [37, E8j for the local-in-time well-posedness results and |40| for linearized stability results of 
stationary solutions. 

The rest of the paper is organized as follows. In Section [2] we reduce the free boundary problem 
(11.1.21) to a fixed domain by using the Lagrangian formulation and state the global existence and 
nonlinear asymptotic stability theorems. The a priori estimates are derived in Section El whose 
main ideas and steps are outlined in Section T3.21 The global existence theorem is proved in Section 
El Section [4] is devoted to the proof of the nonlinear asymptotic stability theorem in the original 
Eulerian coordinates which is an easy consequence of the estimates obtained in Section El The 
local existence of strong solutions in our function space framework given in Section [2] is proved in 
Appendix, Part I. In Appendix, Part II, a linearized analysis is given to illustrate some ideas for 
the original nonlinear problem. 

2 Lagrangian formulation and main results 

2.1 Lagrangian formulation 

In this subsection, we adopt the the Lagrangian particle trajectory formulation as first used in [3j 
for inviscid flows and latter on in [3] [15, [T7] to reduce the original free boundary problem (II.1.21) 
to an initial boundary value problem on the fixed domain x £ [0, R]. 

For this purpose, we first recall some properties of Lane-Emden solutions. For 7 £ (4/3,2), it 
is known that for any given finite positive total mass, there exists a unique solution to equation 
(11.1.51) whose support is compact (cf. [22]). Without abusing notations, x will denote the distance 
from the origin for the Lane-Emden solution. Therefore, for any M £ (0, 00), there exists a unique 
function p(x) such that 

rR 

po := p( 0) > 0, p(x ) > 0 for x £ (0, R) , p (R) =0, M= 4irp(s)s 2 ds-, (2.1.1) 

Jo 

—00 < p x < 0 for x £ (0, R) and p(x) < po for x £ (0, R) ; 

rx 

(p J ) x = —xcj)p, where 4> := x ~ 3 / 4 - 7 T p(s)s 2 ds £ [M/R 3 , 47rpo/3] ; (2.1.2) 

Jo 

for a certain finite positive constant R (indeed, R is determined by M and 7 ). Note that 

(P 7-1 )^ = 7 _1 (7 - 1 ) P ~ 1 ( f) x = -7 _1 (7 - 1)®<£- 
It then follows from (|2.1.1D and (12.1.21) that p satisfies the physical vacuum condition, that is, 

p 7 _ 1 (a;) ~ R — x as x close to R. 
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More precisely, there exists a constant C depending on M and 7 such that 


C 1 (i? — x) < p 1 x (x) < C (R — x) , (0 , R) ■ 


(2.1.3) 


The particle trajectory Lagrangian formulation for (11.1.21) is given as follows. Let x be the 
reference variable and define the Lagrangian variable r{x, t ) by 


rt(x,t) = u(r(x,t),t) for f>0 and r(x, 0) = ro(x), x£l:=(0,R). 
Here ro(x) is the initial position which maps I —>• [0,i?o] satisfying 

(•r o(x) 


rro{x) rx 

/ po(s)s 2 ds = / p(s)s 2 ds, x £ I, 

Jo Jo 

so that 

p 0 (r 0 (x))rl(x)r' 0 (x ) = p(x)x 2 , x el. 


(2.1.4) 


(2.1.5) 


( 2 . 1 . 6 ) 


(Indeed, (12.1.51) means that the initial mass in the ball with the radius ro(x) is the same as that 
of the Lane-Emden solution in the ball with the radius x. Then smoothness of ro(x) at x = R is 
equivalent to that the initial density po has the same behavior near Rq as that of p near R.) The 
choice of ro can be described by 


r 0 (x) = ip 1 (£(x)), 0 <x<R] 
where £ and if) are one-to-one mappings, defined by 

£ : (0, R) —> (0, M) : /' s 2 p(s)ds and : (0, Rq) —> (0, M) : z >->• 

Jo 


(2.1.7) 


/ s 2 p 0 (s)ds. 
Jo 


Moreover tq(x) is an increasing function and the initial total mass has to be the same as that for 
p, that is, 


rRo rro(R) 

/ Airp 0 (s)s 2 ds = / Attpo(s) 
Jo Jo 


rR 

s z ds = / 4irp(s)s 2 ds = M, 

Jo 


( 2 . 1 . 8 ) 


to ensure that ro is a diffeomorphism from I to [0,7?o]- In view of (|1.1.2al) . we see 


l 


r(x,t) 


rr 0 (x) 

p(s,t)s 2 ds = / po(s)s 2 ds, x e I. (2.1.9) 

Jo 

Define the Lagrangian density and velocity respectively by 

f(x,t)=p(r(x,t),t) and v(x,t) = u(r(x,t),t). 

Then the Lagrangian version of (ll.l.2al) and (11.1.2bl) can be written on the reference domain / as 

(r 2 f) t + r 2 f^ = 0 in I x (0, T], (2.1.10a) 

p, ({r 2 v) 3 


(tr\ „ rr 0 (x) 

fv t + -- - + Anfr~ 2 / p 0 (s) 

r x Jo 


s ds = 


r 2 r 7 


in I x (0, T], 


(2.1.10b) 














Solving (I2.1.10al) gives that 


f(x,t)r 2 (x,t)r x (x,t) = Po(r 0 (x))r%(x)r 0 x (x ), x £ I. 


Therefore, 


/OM) = 


x 2 p(x) 


for x £ I, 


r 2 (x,t)r x (x,t ) 

due to (I2.1.6p . So, (11.1.21) can be written on the reference domain / = (0, R) as 


P 



fx 2 /)\ 7 
\r 2 r x J 



v(0,t) = 0, ®(fl,t) = 0 

(r, v)(a;,0) = (r 0 (x), n 0 (r 0 (x))) 


4vr y 2 p(y)dy = p 



X 


in I x (0,T], (2.1.11a) 

on [0,T], (2.1.11b) 

on I x {t = 0}, (2.1.11c) 


where is the normal stress at the boundary given by 

<8 ;= l Al _ A + X2 (*JL + 2 f] = i Al ^ C-) + A 2 - (r2,,) . 


3 r x \rJ x 


r x r 


( 2 . 1 . 12 ) 


Moreover, it can be derived from (12.1.41) . (12.1.5p and (12.1.llbD that 


r(0,t) = ro(0) + / u(0,s)ds = 0 on [0, T], 

Jo 


(2.1.13) 


Notation. Throughout the rest of this paper, c and C will be used to denote generic positive 
constants which are independent of time t but may depend on 7, Ai, A2, M and the bounds of p 
such as p(0) and p(R/ 2); and we will use the following notations: 



:= II • IIl 2 (/), II • \\lp ■= II • \\lp(I) (p = l.oo), and 


| H‘ 


II ' II H a (I) ( s — 1)2). 


2.2 Strong solutions and functionals 

A strong solution to problem (12.1.111) is defined as follows. 
Definition 2.1 v £ C ([0, T]; if ; 2 oc ([0, R))) n C ([0, T]; W 1 ’ 00 ^)) with 


r(x,t ) 



s)ds 


for 


(x,t) £ I x [0, T\, 


( 2 . 2 . 1 ) 


satisfying the initial condition (I2.1.11cl) is called a strong solution of problem ([2.1.Ill) in [0,T] , if 

1) ci < r x (x,t ) < C2, (x,t) £ I x [0,T], for some positive constants c\ and C 2 ; 

2) P~ l/2 [( r 2 r x )~ 1 (r 2 v) x \ x £ C([0, T]; L 2 (I)) and p^ 1 / 2 (r xx , (r/x) x ) £ C ,1 ([0, T]; L 2 (I)); 

3) p 1 / 2 v£C 1 ([0,T]-,L 2 (I)); 

4) v(0,t) = 0 and *B(R,t) = 0 hold in the sense of W 1,ao -trace and H l -trace, respectively, for 
t £ [0, T]; 

5) (12.1.Hal) holds for (x,t) £ I x [0, T], a.e.. 
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Remark 2.2 Let (r,v) be a strong solution of problem (12.1.IIP defined in Definition \2.1l then it 
holds that for any a G (0, R), 


r G C 1 ([0, T]; W 1 ’ 00 ^)) D C 1 ([0, T];H 2 ([0, a])) , r/x G C 1 ([0, T]; H\l)) , (2.2.2) 

u G C ([0,T]; W 1,0 °(/)) rC ([0,T]-H 2 ([0,a])) , v/x G C ([0,T]; H\l)) , (2.2.3) 

(v/r, v x /r x , ®) G C ([0, T]; R 1 (/)) . (2.2.4) 

XTie arguments for (12.2.2P - (12.2.4P go as follows. It gives from c\ < r x (x,t) < c 2 and r(0,t) = 0 
that c\ < x~ 1 r(x,t) < 02 - As a consequence of the facts v G C ([0, T]; W 1,00 (/)) and v(0,t) = 0, 
one has that v/x G C ([0, T]; L°°(I)); which, together with (|2.2.1P and c\ < r/x < C2, gives r/x G 
C 1 ([0, T]; L°°(/)). Moreover, it follows from (|2.2.1P . fT 1-1 ^ 2 (r xx ,(r/x) x ) G C 1 ([0, T]; L 2 (I)), and 
p(x) > p(a) for x G [0, a] that 

(r xx , (r/x) x ) G C 1 ([0,T];L 2 ([0, a])) and (r ra , (v/x) x ) G C ([0, T]; L 2 ([0, a])) . (2.2.5) 

5o, (12.2.2P and (12.2.3P hold, since \(r/x) x \ < a _1 (|r x | + |r/x|) for x G [a, R]. It follows from v/r = 
(v/x)(x/r), ci < r/x < C2, (12.2.2P and (12.2.3|) that v/r G C ([0, T]; id 1 (/)). Dae to of Definition 
\2.1\ it holds that [v x /r x + 2v/r) x G C([0, T]; L 2 (I)). These imply ( v x /r x ) x G C([0, T]; L 2 (/)). ^o, 
(|2.2.4P holds. 


In what follows, we give some remarks on the above definition of strong solutions to problem 
(12.1.IIP . First of all, in order to make the transformation x 1 —> r invertible to define particle 
trajectories, it is essential to require the positive lower and upper bounds for r x , i.e., 1) of Definition 
E31 The requirement v G C([0, +00); kF 1,0O (/)) is also essential, because it is also related to the well¬ 
definiteness of particle trajectories defined by (I2.1.4P . whose uniqueness for the given initial values 
ro(x) is not ensured without the Lipschitz continuity of u in r, and the bound of u r and v x are related 
by the identity u r (r(x,t),t ) = v x (x,t)/r x (x,t). The above definition of strong solutions guarantees 
that each term in equation (|2.1.11a,P is in C([0, T]; L 2 (/)). More than this, each term multiplied 
by p -1 / 2 is in C([0,T];L 2 (/)). Indeed, it is quite natural to require p^ 2 v t G C([0, T]; L 2 (I)) other 
than pvt G C([0, T]; L 2 {I )) from the kinetic energy point of view. The kinetic energy of the system 
is ||xp 1 / 2 u(-,f)||, which is expected to be bounded and continuous in time. (Indeed, this can be 
justified by applying the multiplier r 2 v to equation (I2.1.11ap .) By studying the problem obtained 
by differentiating the original one with respect to t, one may expect ||xp 1 / 2 ut(-, t)\\ is bounded and 
continuous in time. We can improve this to the boundedness and continuity in time of ||p 1 / 2 wt(-,f)|| 
by using the viscosity, which will be shown later. This leads us naturally to require that each term 

in (|2.1.11a,P multiplied by p -1 / 2 is in C([0, T]; L 2 (I)). _ 

To show the well-posedness of strong solutions to problem (12.1.IIP defined in Definition 12.11 we 
introduce the following higher-order functional: 


£(f) 


1) v x)('i t) ||o + 


fA~ 1 / 2 (r x 


(r/x) x ){-,t) 


+ 


P 1/2 vt(-,t ) 


2 , t> 0. (2.2.6) 


We will prove the global existence and uniqueness of strong solutions satisfying <£(t) < C(£(0) 
(t > 0) for some constant C > 0 independent of t, and some decay estimates as in i) of Theorem 
13.11 provided that (£(0) is suitably small and the following compatibility condition of the initial 
data with the boundary conditions holds. 


u(0,0) = 0 and Q3(i?, 0) = 0. 


(2.2.7) 


10 










































In order to gain further regularity of strong solutions obtained via <£(f), we introduce the 
following functional: 


$a{t) 


«(*) + 


_(2 7 -l-a)/2 ? 




a £ [0, 2y — 1], t > 0. 


Besides £(0) is small, if the initial data are assumed to satisfy further regularity that 3a(0) < 00 
for 0 < a < 2y — 1, we will prove further regularity and decay (than those stated in 1 ) of Theorem 
I3TD of strong solutions with decay rates which may depend on a in various norms, as shown in ii ) 
and in) of Theorem O Indeed, the index a indicates the behavior of solutions near the vacuum 
boundary, which has influence on the decay rates of solutions. In particular, we will prove the 
regularity that v xx (-,t) £ L 2 (J) for all t > 0 and the decay of ||u X x(')i)|[> if the initial data satisfy 
additional regularity that r xx (- ,0) £ L 2 (I) (i.e., 32 7 -i(0) < 00). 

Some remarks are given on the smallness requirement of functional £(t). For problem (12.1.111) . 
r(x,t ) = x is the equilibrium solution. The basic stability requirement is on the smallness of 
r(x,t ) — x for all time in certain topology. The smallness of the quantity r x (x,t ) — 1, the x- 
derivative of r(x,t) — x , in L°°-norm, ensures that r x (x,t) is bounded from below and above by 
positive constants. This also gives the smallness of \\r/x — l||i,°° and thus the lower and positive 
bounds for r/x due to the condition r(0,f) = 0 (see (12.1.131) 1. So, it is essential to obtain the 
smallness of \\r x — 1||l°°. Since v = 0 is the equilibrium for problem (12.1. Ill) , it is required that 
||r>(-, Ollw 1 ’ 00 ^/) f° be small for all t > 0 if this holds true initially, from the stability point of view. 
It is crucial to derive the smallness of T°°-bound for v x , which also gives the bound for v/x due to 
the condition that v (0,f) = 0. Our basic strategy is to use the weighted T 2 -energy method together 
with some pointwise estimates to bound ||r x — 1 ||l°° and H^H^oo. For example, we may bound 
| \r x - 1||l°°([o,.r/ 21) use of \\ r xx\\L 2 {[o,R/ 2 ]) which can be bounded by ||p 7 ~ 1/2 r a;a; ||. Therefore, that 
the functional (£(t) is small for all time is essential for our stability analysis of problem (12.1. 111) . 

We give here some consequences of the smallness of (£(t) (f > 0), which ensure that each term 
in (12.1. Hall multiplied by p^ 1 / 2 is in L 2 (I ) (t > 0), and that the boundary condition v(0,t) = 0 
and i8(R,t) = 0 are achieved in the sense of IF 1,00 -trace and i7 1 -trace, respectively. For t > 0 and 
a £ (0, R ), it holds that 


II (r/x - l)(-,t)||ioo < ||(r* - l)(-,t)||i<» < £(i), 

1/2 < r(x,t)/x < 3/2 and 1/2 < r x (x,t ) < 3/2, x £ I, 

\\{v/x)(-,t)f Loo < |M-,i)||ioo < C(t), 

P~ l/2 (v x /r x + 2v/r) x (-,t) = p~ 1/2 [{r 2 r x )~ 1 (r 2 v) x ] x (-,t) <C£(t), 

P y ~ 1/2 (v xx , (v/x) x )(-,t) < C£(t), 


(v,r - x)(-,t)\\ 2 H2{[0 a]) < C(a)<£(t), 


-- 1,-,-A® ) (*,*) 

x x r r x 


H 1 


< cm, 


( 2 . 2 . 8 ) 

(2.2.9) 

( 2 . 2 . 10 ) 

( 2 . 2 . 11 ) 

( 2 . 2 . 12 ) 

(2.2.13) 


where C and C(a ) are positive constants independent of t. 

The arguments for ()2.2.8l) - (12.2.13j) go as follows. (|2.2.8I) and (12.2.101) follow from r(0,t) = 0 and 
v(0,t) = 0, respectively; (12.2.91) follows from the smallness of <£(t); (12.2.111) follows from (12.2.81) . 
(12.2.91) . the definition of £(t) in (12.2.61) . and equation (I2.1.11al) which can be rewritten as 




(' r 2 v) x 

r 2 r x 


= P (r) Vt + p1 


x 2 1 
r 2 r x 


+ 


x 2 1 \ 7 x 4 " 1 

r 2 r x ) r 4 


(m 
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(12.2.121) follows from (12.2.101) . (12.2.111) . the definition of <£(t) in (12.2.61) . and the estimate below: 



(Indeed, (12.2.14)) follows from (I3.4.5|) . (13.4. 7|) . (12.2.91) . and (|3.4.23j) which will be prove later.) In a 
similar argument as that for Remark 12.21 one can obtain ()2.2.13p . 


2.3 Main theorems and remarks 


The first theorem of this paper is on the global existence and the regularity of strong solutions, 
which also gives the strong Lyapunov stability of the Lane-Emden solution in the functional 


Theorem 2.3 Let 7 £ (4/3, 2) and p be the Lane-Emden solution satisfying (|2.1.1D - (|2.1.2D . As¬ 
sume that (12. 2. 71) holds and the initial density po satisfies (11.1.31) and (|2.1.8D . There exists a 
constant 5 > 0 such that if (£(0) < 5, then the problem (12.1.Ill) admits a unique strong solution in 
I x [0, 00) with 

<S(t) < C<£(0), t > 0, (2.3.1) 


for some constant C independent oft. 

Moreover, i/52 7 -i(0) < 00 (i.e., ||r xx (-,0)|| < 00 ), then the strong solution obtained above 
satisfies the following further regularity estimates 

Kx(-,i)|| 2 <C5 27 -i(0)(1 + 5 27 -i(0))(1 + C(0)), t> 0, (2.3.2) 

for some constant C independent of t; and 

\\r xx (;t)\\ 2 < C52 7 —i(0) + C(T)(£(0), t £ [0,T], (2.3.3) 

for some constants C independent oft and C(T ) depending on T. 

For any t > 0, since r x (x,t ) > 0 for x £ L. r(x,t ) defines a diffeomorphism from the reference 
domain I to the changing domain {0 < r < R(t)} with the boundary 

R(t) = r(R,t). (2.3.4) 

It also induces a diffeomorphism from the initial domain, B Rq ( 0), to the evolving domain, -£>777 (0), 
for all t > 0: 

x + 0 £ B Ro ( 0) -»• r (r„ 1 (|x|),t) |x| _1 x £ B R{t) ( 0), 
where rfi 1 is the inverse map of vq defined in (|2.1.7I) . Here 

B r o(0) := {x £ M 3 : |x| < Rq} and B R ( t )( 0) := {x £ M 3 : |x| < R(t)}. 


Denote the inverse of the map r{x, t) by Lit f° r t > 0 so that 

if r = r(x,t) for 0 <r<R(t), then x = lZt(r). 

For the strong solution (r,v) obtained in Theorem 12.31 we set for 0 < r < R(t) and t > 0, 
x 2 p(x ) 


p(r,t) = 


r 2 (x,t)r x (x,t) 


and u(r,t) = v(x,t) with x = Rt(r). 


(2.3.5) 


Then the triple (p(r,t),u(r,t), R(t)) (t > 0) defines a global strong solution to the free boundary 
problem (11.1.21) . Furthermore, we have the strong nonlinear asymptotic stability of the Lane-Emden 
solution as follows. 
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Theorem 2.4 Under the assumptions in Theorem \2Jft Then the triple ( p,u,R(t )) defined by 
(|2.3.4I) and (|2.3.5j) is the unique global strong solution to the free boundary problem (11.1.11) satisfying 
R G JT 1 ,<x, ([0, +oo)). Moreover, the solution satisfies the following estimates. 

i) For any 0 < 0 < 2(7 — l)/(37), there exists a positive constant C{9) independent of t such 
that for all t > 0, 


sup \r(x,t) — x\ < C(9)(l + t) h >VW), 

0 <x<R 

sup \u(r,t)\ < C(6)(l + t)~ + 2y / G:(0), 

0 <r<R(t) 

sup \{u r , r~ L u) (r,t)\ < C(9)(l + t) W + 2 

0 <r<R(t) 


sup 

0<x<R 


(p(x)f 7 ~^ 4 [p(r{x, t),t) - p(x)] < C(9)( 1 + ty 


7-1 
2 7 " 


: 7e(o). 


(2.3.6) 

(2.3.7) 

(2.3.8) 

(2.3.9) 


ii) Suppose that 3a(0) < 00 for some a G (0,7). Let 6 be any constant satisfying 0 < 9 < 
min {2(7 — 1)/(37), 2(7 — ct)/7}. Set k = (I/ 7 ) min{o — (7 — 1), 7 — 1} — 9 when a G (7 — 1, 7 ), 
and k = 0 when 0 = 7 — 1. Then there exists a positive constant C{a , 9) independent oft such that 
for all t > 0, 


sup 

0<x<R 


(p(x)) {37 6 a)/4 [p(r(x,t),t) - p(x)\ < C(a,9)(l + t) ^ +2 \AM0); 


and if a G [7 — 1, 7 ), 

sup \u(r,t)\<C(a,0)(l + t) ^ ^/Sa(0) + <£(0)5 a (0), 

0 <r<R(t) 

sup | (u r , r~ 1 u)(r,t)\ < C(a,9){\ + + <£(0)&,(0), 

0 <r<R{t) 


sup 

0 <x<R 


{p{x)) (l 2)/2 {p(r(x,t),t) - p{x)} 


< C(a,9)(l + t )-f-^ 1 +f ^(0) + C(0)ff Q (0). 


Here 


b-i = min < max 


2 27 2 


a + 1 3 27-I 9 

-. —R + - — — 

27 -lla’ 2 27 2 


k 47 — 3 3 

- + —- 1 

2 27 2 


U-lM-e, 

7 


k 47 — 3 3 r IO 7 — 9 9 


bo = min < —h -, , 

4 2 27 2 ’ 4 4 7 


+- + 


k 47 — 3 3 


2 7 


(2.3.10) 

(2.3.11) 

(2.3.12) 

(2.3.13) 


(2.3.14) 

(2.3.15) 


Hi) Furthermore, if \\xp l / 2 vtt(-, 0) || 2 + \vt(R, 0)| < oo, then R G TT 2 ,oo ([0,+ 00 )) and 
TOI < M^,0)| +C(^(0 )) 1 / 4 ((e(O )) 1 / 4 + \\ X p l / 2 v u (;ti)\fi) , t > 0, (2.3.16) 

for some constant C independent oft. 
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Remark 2.5 In <£(0), the term \\p 1 / 2 vt(-,0)\\ can be given in terms of the initial data as follows: 
for x G I, 


P*Vt(x, 0) = P a Q |At ^ 


(r 2 v) a 


r * r 


-P 1 


x 2 1 


r 7 


x 2 1 


r- r~ 


x 


(. P y )x ( (®>0)- 


Indeed, this is equivalent to equation (I2.1.11al) at t = 0. 


Remark 2.6 The estimates in (|2.3.9|) . (|2.3.10|) and (12.3.131) yield the uniform convergence with 
rates of the density to (11.1.11) to that of the Lane-Emden solution for both large time and near the 
vacuum boundary since 7 < 2. 


Remark 2.7 The initial perturbation here includes three parts: the deviation of the initial domain 
from that of the Lane-Emden solution, the difference of initial density from that of the Lane-Emden 
solution, and the velocity. Since the Lane-Emden solution is completely determined by the total 
mass M, our nonlinear asymptotic stability result shows that the time asymptotic state of the free 
boundary problem is determined by the total mass which is conserved in the time evolution. 

Remark 2.8 The condition \\xp l / 2 vu{-, 0) || 2 + \vflR, 0)| < 00 in in) of Theorem \2.f\ to ensure 
R G W 2,oo ([0, +00)) (uniform boundedness of the acceleration of the vacuum boundary) is a higher- 
order compatibility condition of the initial data with the vacuum boundary. Indeed, one may check 
from the proof that every particle moving with the fluid has the bounded acceleration for t G [0, 00) 
if it does so initially. 


3 Proof of main results 


3.1 A theorem with detailed estimates 

For the convenience of presentation, we set R = 1 and I = (0, R) = (0,1). Indeed, we will prove 
the following results for the global strong solutions obtained in Theorem 12.31 which gives not only 
the nonlinear asymptotic stability results stated in Theorem 12.41 but also detailed behavior of the 
solutions both in large time and near the vacuum boundary and the origin. 


Theorem 3.1 Let v be the global strong solution to the problem ((2.1.111) with r given by (12.2.lj) as 
obtained in Theorem \2.3l 

i) Let 9 and 5 be any constants satisfying 0 < 6 < 2(7 — l)/(37) and 6 G (0,1). Then there 
exist positive constants C{9) and C{9,5) independent oft such that for all t > 0, 


2 ( 7 — 1 ) 


(1 + t) 

+ (1 +t) 


(r - x)(-,t)\\ 2 Loo + (l + t) 27 \\(v,xv x )(-,t)\\ 2 Lo 


A 0 p( 3 7 2 )/ 4 (r x - !,r/x _!)(•, t) +(1 + A) : V 6 


( v x ,v/x ) (-,t)||io 


+ (1 + t) y 6 ^ (xp 1/2 v t ,v,xv x ^ (-,f) + p 7/2 (r - x,xr x - x) (-,f) ^ 


3(7-1) 


+ (l + t) 7 

+ 


-e 


r-x)(-,f)|| 2 + (l + t)V 6 (r x -l,r/x-l,v x ,v/x,p 1/2 v t ^{-,t) 


p- * 7 2 (r — x, xr x — x) (•, t) < C(0)£( 0) 


and 


(1 + t) 1 r L - e 


(r x - 1 ,r/x - l,v x ,v/x) (•, 11 ifi ([o,<5]) < C(9,5)<£( 0 ). 


(3.1.1) 

(3.1.2) 
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ii) Suppose that 3a (0) < oo for some a £ (0,7). Let 6 and 6 be any constants satisfying 
0 < 6 < min{2(7 — l)/(37), 2(7 — a)/and 5 £ (0,1). Then there exist positive constants 
C{a), C(a,9) and C (a, 9,5) such that for all t > 0. 


da(t) < C(a)da( 0 ), 


(1 +f)^- 1)h ~ e 
and if a £ [7—1,7), 


37 2 a)/A {r x - l,r/x - l)(-,t) <C(a,6)$ a (0y, 

L°° 


(3.1.3) 

(3.1.4) 


K I 47—3 3 n 

(1 + t ) 2 2 T 2 


v x ,v/x,p 1/2 v t J (-,t) +(l + t) ^ +4 4 ° ||v(-,t)||£<, 

+ (1 +t)fr\\xv x (.,t)\\l- + (1 +i)^ min{bl ’ b2} \\(v X ,v/x) (;t)\\l~ 


87 — 5 1 k _ 5n 

\ Arv I A A U 


(3.1.5) 


+ (1 +£) 2+ 2 t 2 f /2 {r x - l,r/x - l)(-,i) <C(a,0)5 a (O), 

L°° 

min f K I 47-3 3 C it I 107-9 9q1 „ 

(1 + t) 1 2 27 2,4 47 4 S \\(r x -l,r/x-l,v x ,v/x)(-,t)\\ H i 


([0,5]) 


K I 47—3 3 Q 


(3.1.6) 


+ (l + t) 2 27 2 11(77 - l,r/x - 1 ) (-,t)|| L 2 (M) < C(a,6,5)$ a (0). 


Here 3 a (0) = 3a(0) + £(0)3a(0), k = 0 when a = 7 — 1, k = (I/7) min{o; — (7 — 1 ), 7 — 1} — 6 
when a £ (7 — 1,7), and b\ and 62 are given by ()2.3.14l) and (12.3.15|) . 

in) Suppose that 3a(0) < 00 for some a £ [ 7 , 27 — 1], Let 6 be any constant satisfying 0 < 6 < 
2(7 — l)/( 37 ). Then there exist positive constants C and C(9) such that for all t > 0, 


3a(t) < Cda( 0) + C(6)( 1 + f)(“-' r + 0 T , )/( a - 1 )g( O ); 

and if a = 27 — 1. 

Ik^(-,t)|| 2 < ^27-1(0) + <7(0X1 + t)5 + 3^g( 0 ), 


77-6 


V XX , (v/x) x ) < C(0)( 1 + t)-^r+ 4e ^_ 1(0) (1 + fo 7 -i(0)) (1 + 0(0)). 


(3.1.7) 

(3.1.8) 

(3.1.9) 


iv) Suppose that ||cc / o 1//2 ^it(-, 0)|| 2 < 00. Then there exists a positive constant C independent of 
t such that for all t > 0, 

poo 

\\xp 1/2 Vtt{-,t)\\ 2 + / \\{v ss ,xv ssx )(-,s)\\ 2 ds < C£(0) + C\\xp 1/2 v u (-,0)\\ 2 . (3.1.10) 

J 0 


3.2 Main ideas and the structure of the proof 


The local existence and uniqueness of strong solutions to (12.1.111) on some time interval [0,T*] are 
given in Appendix, Part I, by using a finite difference method as used in [321 25, :26j [2j. In order to 
prove the global existence of strong solutions, we need to derive the uniform-in-time boundedness of 
the nonlinear functional (£(i) defined in ([2.2. 61) . Our basic strategy is to use the weighted L 2 -energy 
method together with some pointwise estimates. For the weighted L 2 -energy estimates, motivated 
by the linearized analysis for problem (12.1.lip shown in Appendix, Part II, a natural functional 
should be £{t) given by 


£(t) := || (r - x,xr x - x)(-,t)|| 2 + \\{v,xv x ) (-,t )|| 2 + IKr* - 1 )M)||£<» ([1/2i 1]} 


+ 


P 7 l/2 (r xx , (r/x) x )(-,t) + p 1/2 v t (-,t) 


(3.2.1) 
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(Indeed, (£(t) and £{t) are equivalent under some assumptions which will be shown in Lemma [3. 171 1 
However, much efforts are needed to pass from the linear to nonlinear analysis. An important step 
for this is to identify the a priori bounds with which the basic bootstrap argument can work. It is 
found that the appropriate a priori assumption is that | r x (x,t) — 1| and |u x | are suitably small to 
pass from the linear to nonlinear analysis. To this end, we use a bootstrap argument by making 
the following a priori assumptions. Let v be a strong solution to (12.1.111) on [0, T] with 

r(x,t ) = r 0 (x) + f v(x,t)cIt, (x,t) G [0, 1] x [0, T}. 

Jo 

The basic a priori assumption is that there exist suitably small fixed constants eo G (0,1/2] and 
e\ G (0,1] such that 

\r x (x, t) — 1| < eo and \v x (x, t)| < e\ for (x, t) G I x [0, T]. (3.2.2) 

It follows from (13.2.21) and the boundary condition u(0, t) = 0 (so r(0, t) = 0) that 

| x~ 1 r(x,t) — l| < | r x (x,t) — 1| < eo for (x,t) G [0, 1] x [0,T], (3.2.3) 


|x 1 v(x,t)\ < \v x (x,t)\ < e\ for (x,t) G [0, 1] x [0, T], (3.2.4) 

In particular, it holds that 

1/2 < r x (x,t ) < 3/2 and 1/2 < x~ 1 r(x,t) < 3/2 for (x,t) G I x [0, T], (3.2.5) 


In order to close the argument, it is needed to bound the L°°-norms of r x — 1 and v x by the 
initial data. The usual approach for this is to use energy estimates and the Sobolev embedding, for 
example, for r x — 1, 


IIO* - l)(-,«)l|!» < c (||(r, - 1)(-,()|| 2 + ||r„(-,t)|| 2 ) . 


However, due to the strong degeneracy of the equation near the vacuum boundary, the L 2 -norm of 
r xx may grow in time (see (13.1.81) ). the uniform L 2 -bound for r xx is valid only for an interval of x 
away from the vacuum boundary x = 1 (see (|3.1.2I) L say, Hycca;||z, 2 ([o,i/2])- The term \\p^^ l / 2 r xx {-,t)\\ 
in the functional £{t) can only give the bound of | r x (x, t) — 1| for x away from the vacuum boundary. 
This is the reason why we include ||(r*^ — 1)(*,i)([i /2 i]) ' n ^ ie functional £(t). Similar ideas 
apply to the estimate of ||ua;(-, t )||loo , which is routinely bounded by ||u x (-, t )||//i via energy methods 
in the region away from the vacuum boundary x = 1, say [0,1/2], since the system is not degenerate 
there. However, in the region near the vacuum boundary, say, [1/2,1], it is not a easy task to obtain 
the L 2 -bound of the second derivative of v, ||uxx(-,t)||i2/r 1 / 2) i]), again, due to the strong degeneracy 
of the system near the vacuum. Indeed, our idea to bound v x near the vacuum boundary is to 
use the L 2 -norm of the viscosity term to obtain the L 2 -bound for (v x /r x ) x , instead of that for v xx . 
Because r x has positive lower and upper bounds. 

The L 2 -norm of the viscosity term is quite different from that of v xx . This can be seen as 
follows: 


V : = 


{r 2 v) x 

r 2 r x 


X 



9 v x r xx + 2 ( 
r x V r / x 


Under the a priori assumption that r x is close to 1, the difference between the viscosity term V 
and v xx involves the second derivative of r, r xx , whose L 2 -norm may grow in time. This growth 
may be balanced by the decay of the L°°-norm of v x so that the L 2 -norm of r x 2 v x r xx is bounded. 
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It is possible to make an ansatz of the decay of the L°°-norm of v x and then close the argument. 
However it is quite intricate to do so. On the level of the second derivatives, our strategy is to close 
the estimates by using the L 2 -estimate of the viscosity term V, instead of that of v xx . It should be 
emphasized that it is enough to get the L 2 -estimate for the viscosity term V to close the estimates 
and obtain the necessary bounds for the global existence of strong solutions and assurance of the 
well-definiteness of the boundary condition, as explained in Section 12.21 After proving the global 
existence of strong solutions and basic decay estimates, we can improve further regularity that 
v xx G L 2 (I) for all t > 0 and the decay of ||v X x(-,i)|| if the initial data satisfy additional regularity 
that r xx (-, 0) G L 2 (I) (i.e., £27-1(0) < 00). 

Due to equation (13.2.711 , the equivalent form of equation (12.1.Hall , the L 2 -norm of the viscosity 
term V can be bounded by C£(t). (More precisely, ||/3~ 1 / 2 V(-,f)|| < C£(t).) So, it suffices to show 
that the higher-order functional £(t) defined by (|3.2.ip is bounded uniformly by the initial data, 
i.e., 

£(t)<C£( 0) for all t G [0.T], 

We outline the main steps for this as follows. 

Lower-order estimates. The key elements in our analysis are the weighted estimates by applying 
various multipliers to the following equation: 


P 



f x 2 p 
\r 2 r x 



®x+4A;l(-) , 
\r J x 


which is equivalent to (12.1. llal) . Here IB is defined in (12.1.1211 . In Lemma l3.5l we use the multiplier 
r 2 v to yield the bound for the basic energy 


xp 1/2 v(-,t ) + xp l/2 (r/x - l,r x - 1) (-,t) +( ||(u, xv x )(-, s)|| 2 ds. 

Jo 


The multiplier r 3 — x 3 plays an important role in the proof of Lemma 13.61 in which a bound is 
obtained for 


x (r/x - 1 ,r x 



i)M 


2 

ds, 


which refines the weighted estimate of \\xp^^ 2 (r/x — 1 ,r x — 1) || obtained in the basic energy es¬ 
timates. We also show the decay estimates for the basic energy in Lemma 13.61 by establishing a 
bound for 


(1 + 1 ) 


xp l / 2 v(-, t) 


2 

+ 



_ np rp rp 

Jjl x 


x) (-,t) 


+ [ (1 + s) ||(u,xu x )(-,s)|| 2 ds. 
Jo 


With those estimates, we are able to bound \r x — 1| away from the origin in Lemma 13.71 by noting 
the fact that 53 can be written as the time derivative of a function so that one can integrate equation 
(13.3.111 with respect to both x and t to get the desired estimates, where the monotonicity of the 
Lane-Emden density plays an important role. Adopting the multiplier r 2 vt, a bound for 


(1 + 1 ) 


xp 1/2 v t (-,t ) 


+ \\(v,xv x ) (-,t)|| J ^ “I" J 0 (1 + s ) IK u *’ x Wz)(-,s)|| 2 ds 


is given in Lemma l3.8l bv studying the time differentiated problem of (12.1.1111 . 

Further decay estimates and regularity are given in Lemma 13.91 which is important to the 
derivation of the decay for ||r — x||loo in (I3.1.1H . This in particular implies the convergence of the 
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evolving boundary r = R(t ) to that of the Lane-Emden stationary solution. Lemma [3.91 also shows 
that rates of time decay in various norms depend on 7 indicating the balance between the pressure 
and self-gravitation, and the dissipation of the viscosity. These further decay estimates are derived 
from the following two multipliers: 

px px 

/ P~^(y)(r 3 - y 3 )ydy and / p~ l3 (y)(r 2 v)ydy for 0 < ft <j -1. (3.2.6) 

Jo Jo 

It should be noted that the first multiplier in (13.2.61) is motivated by virial equations in the study 
of stellar dynamics and equilibria (cf. [2(71 [39]). To the best of our knowledge, those multipliers 
have not been used in previous literatures. 

Higher-order estimates. To get the L 2 -estimate of the viscosity term V, we write it as follows 


V = Q xt , where Q := In r x + 2 In - , 


and rewrite (I2.1.11ap as, by virtue of (12.1.21) . 

X 2 \ 7 / x \ 4 

-pvt- ( - 5 — ) - (- 

r 


pGxt +7 


x 2 pV _ x 2 _ 

9 ) z/x — 9 P^t 

r z r I 


r 2 r. 


xcf>p , where (j>(x) = x 6 Aitp(s)s z ds. (3.2.7) 

Jo 


This form has advantage in estimating the second derivative (in x) estimates here. Indeed, the 
interplay among the viscosity, pressure and gravitational force can be seen easily. The gradient of 
the pressure is decomposed as follows: 




7 

xcpp — 7 



Qx- 


The first part of this decomposition is used to balance the gravitational force and the second part 
is the t- antiderivative of the viscosity, Q x t, multiplied by a weight which is equivalent to p 7 . This 
weight is degenerate on the boundary, but strictly positive in the interior. The degeneracy near 
the vacuum boundary is one of main obstacles in higher-order estimates, which is overcome by 
choosing suitable weights and multipliers as we will outline below, and a delicate use of the Hardy 
and weighted Sobolev inequalities. (The method of using Hardy and weighted Sobolev inequalities 
to build regularity was first adopted for the physical vacuum problem for inviscid flows in [3], and 
latter on in Husnu].) Indeed, in terms of Q, the principal part of (13.2.71) is 

/ x 2 \ 7 

pGxt + 7 ( -9— ) P^Gx, 

\r-r x J 


which is linear in Q x and with a degenerate damping. This structure leads to desirable estimates 
on Q and their derivatives. 

With the lower-order estimates obtained already, we can derive in Lemma 13.181 the uniform 
bound for ||p 7_1 / 2 (r a;a ;, (r/x) x ){-, t)|| (due to the bound for ||p 7 ~ 1 ^ 2 I7x( , T)|| and (13.4.22]) ) and the 
decay estimates for 11p 1 / 2 ^^(-, t) ||. This completes the proof of the uniform-in-time bounds for the 
higher-order energy functional £ ( t ), which also verifies the a priori assumptions (13.2.31) and (13. 2. 41) 
due to the equivalence of £(t) and <£(t) shown in Lemma l3.17l and consequently, the global existence 
of the strong solution is obtained. 

With the decay estimates for the lower-order norms in Lemma 1,3.91 and the higher-order estimates 
in Lemma [3.181 we prove the decay estimates of ||(r x —1, r/x — 1) (-, t)||, ||(u, v x )(-,t) || L <x> and 
||(r — a:)(-, t)\\ Loo in Lemma 13.191 with which part i) of Theorem 13.II is proved. 
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The second part of the higher-order estimates will be given in Section T3. 4. 31 in which the faster 
decay estimates are given under the assumption of the finiteness of 3a(0), a G (0,7). A key 
ingredient in the proof is to use the new multiplier x 2 p 2 ' y ~ 2 ‘y$ti where 


= 7 


x 2 p 

r 2 r T 


Qx + 


r 2 r. 


x 4 > Pi 


the sum of the gradient of the pressure and gravitation force. The proof of ii ) of Theorem 13.11 
consists of Lemmas I3.2UH3.221 

Parts in ) and iv ) of Theorem [3TT1 on the further regularity of solutions are proved in Section 
13.4.41 In particular, the bound for \\r xx (-, t)|| and the decay for ||u xa; (-, t)\\ are given in Lemma [3.241 


Remark 3.2 The insight for the above estimates and the basic decay rates can be gained from 
the linearized analysis for problem (12.1.111) . which is shown in Appendix, Part II. The linearized 
analysis gives the ideas for the basic decay rates. The further decay rates of various norms in 
Theorem \3.1\ (based on Lemma 1 3 . 511 are obtained by applying various multipliers, reflecting the 
fact that the asymptotic stability mechanism is due to the balance between the pressure and self¬ 
gravitation, and the dissipation of the viscosity, and also the decay rates depend on the behavior of 
the initial data near the vacuum boundary. In particular, for the multipliers in (13.2.61) . we choose 
/3 < 7 — 1 so that f p-^dx < 00 due to the physical vacuum behavior that p{x) ~ (1 — x) 1 /^ 7-1 ) 
near the vacuum boundary. This has an effect on the decay rates of various norms given in Lemma 
14. .91 which depend on 7 , and indicates that the physical vacuum behavior may be the reason for the 
dependence of decay rates on 7 . 


In the rest of this article, we will frequently use the following weighted Sobolev embedding and 
general version of the Hardy inequality, whose proof can be found in [8], to build the regularity. 
This idea was first used in [3J for inviscid flows, and latter on in smani]. 

Lemma 3.3 (weighted Sobolev embedding) Let d denote the distance function to the boundary dl. 
Then the weighted Sobolev space H d (I), given by 

H d (I) := I dFe L 2 (I) : d 2 (|T| 2 + |F X | 2 ) dx < 00 j , 

satisfies the embedding: 

H l d {I)^L 2 {I). (3.2.8) 


Lemma 3.4 (Hardy inequality) Let k > 1 be a given real number and g be a function satisfying 

, 1/2 


r-L/Z 

/ x k ( g 2 + g 2 x ) dx < 00 , 
■Jo 


then it holds that 

r-l/2 


/ 


1/2 


x k 2 g 2 dx < c / x k (g 2 + g 2 ) dx, 


(3.2.9) 


where c is a generic constant independent of g. 


As a consequence of Lemma 13.41 one has 

[ (1 - x) k ~ 2 g 2 dx < c [ {l-x) k (g 2 + g 2 x ) dx, 
J 1/2 J 1/2 

provided that the right-hand side is finite. 


(3.2.10) 
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3.3 Lower-order estimates 


In this and next subsections, we derive the a priori estimates for the strong solution (r, v ) on the 
time interval [0 ,T] defined in Definition 12.11 under the assumption (13.2.311 and (13.2.41) . We start 
with the lower-order estimates in this subsection, for which we rewrite equation (12.1.llal) as 


v t + 


2 — \ T 

x z p x 

r 2 V x 


--T + 4Ai - 


7V x 


Here 23 is defined in (12.1.121) . Recall the boundary conditions (12.1.11bl) and (12.1.ID : 
p(l) = 0, v(0,t) = 0 and 23(l,f) = 0. 

First, we estimate the basic energy, for which the condition 7 > 4/3 is crucial. 
Lemma 3.5 Suppose that (13. 2. 3D holds for a suitably small positive number eo- Then, 

xp*v(-,t) + (3y — 4) xp* - l,r x - l) (-,t) +cr f \\{v,xv x ) (-,s)|| 2 ds 


(3.3.1) 


(3.3.2) 


<c 


xpiv(-,0) 


+ 


-4 f r ° 1 1 

Xp 2 - 1 , r 0x ~ 1 

V x 


(3.3.3) 


0 < t < T, 


where a = min{2Ai/3, A2}. 


Proof. Multiplying equation (I3.3.1|) by r 2 v and integrating the product with respect to spatial 
variable, we have, using the integration by parts and boundary condition (13. 3. 2D . that 


j fj(x,t)dx = — J 23 ( r2 v) x dx + 4Ai j r 2 v (-) dx, 


where 


fj(x,t ) := ^x 2 pv 2 +x 2 p ) 


1 fx \ 2 7-2 /1\ 7 /®\ 2 .a; 

' ' — + ( -) r x - 4- 

r x J V r / r 


4 


7 — 1 V r 

By the Taylor expansion, the quantity [■] in 77 can be rewritten as 


(3.3.4) 


(3.3.5) 


4-37 

7-1 


+ (2 "^(" r 0 + 2 (^ _1 ) +(r *- 1 ) 


+ Q, 


where Q represents the cubic terms which can be bounded by 

3 n 


|Q| < c | r x - 1| + 


- - 1 

x 


< ce 0 \r x - 1| + 


r 

- - 1 

x 


due to (|3.2.5D and (|3.2.3D . This implies that for 7 G (4/3, 2], 


1 / x \ 27—2 / 1 


7 — 1 V r 


7-1 


x 


— + - r x - 4- > 


x 4 — 37 37 — 4 


7-1 


+ 


2 (/- 1 ) 2 + h - 1 ) ! 


provided that eo is less than a constant depending on 37 — 4. Set 
rj{x,t) := fj(x,t ) - -—^x 2 p 7 . 

7-1 


(3.3.6) 
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Then the above calculations imply that 


1l(x,t) > \x 2 pv 2 + 37 4 


2 G- l ) 2+ ( r «- l ) s 


(3.3.7) 


T](x,t ) < -X 2 pV 2 + CX 2 p 1 


(l -0 a +(--‘) s 


Clearly, (|3.3.4|) and (12.1.121) show that 


d f , , , 4 7 r 4 

it. J „(x,t)<fc = --V j - 


V 

r / x 


dx-X 2 ^ 2 \{r 2 v)fdx. 


(3.3.8) 


(3.3.9) 


Note that 


r 

r x 


v\ 2 r 2 


r/ x 


= —u 2 + r^u 2 — 2rvv x and 


(r 2 n 


■) C „2 


r r r^ 


31 = —v 2 + 4r x v 2 + 4 rvv x . 


We obtain 


d_ 

dt 


J rj(x,t)dx < —3a J 


-v x + 2 r x v 


dx, 


where a = min{2Ai/3, A2}; and 


(3.3.10) 


— v 2 + 2r x v 2 dxds < I rj(x,0)dx, t £ [0,T], 


(3.3.11) 


j rj(x, t)dx + 3a u 
This, together with (|3.3.7D . (13.3.81) and (|3.2.5D . implies (13.3.31) . 

□ 

In the following lemma, we use the multiplier r 3 — x 3 which is motivated by virial equations 
in the study of stellar dynamics and equilibria (cf. [201 39]) to refine the weighted estimate of 
|| xp 1 ! 2 (r/ x — 1, r x — 1) || obtained in Lemma l3.5l bv improving the estimates near the vacuum, and 
give the decay estimates for the basic energy. 

Lemma 3.6 Suppose that (13. 2. 3D holds for a suitably small positive number eo- Then, 


a 


<C 


X (x - 1,Vx ~ 1 ) ('>*) + ( 3 7 — 4 ) / xp^ - l,r x - l) (-,s) 


ds 


/ r ° 1 1 

x( - 1 , r o x — 1 

X 


+ 


xp 2 v(-, 0) 


(3.3.12) 


, 0 <t<T, 


and 


(1 + t) xpiv(-,t) + (3y — 4)(1 + t) xp 2 - l,r x - l) (-,t) 

+ (l+t) j w ||(r-x)(-,t)|| 2 + a / (1 + s) \\{v,xv x ) (-,s)|| 2 ds 

Jo 

<c(^x( 


(3.3.13) 


11 xq x 1 


xp 2 v(-, 0) 


, 0 < t < T, 


where a = min{2Ai/3, A2}. 
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Proof. The proof consists of two steps. With the basic energy estimate obtained in the previous 
lemma, we can achieve the estimate for \\x(r x — l,r/x — 1|| by a moment argument in Step 1. It 
should be pointed out the double integral obtained in Step 1 will play a crucial role in the derivation 
of higher-order estimates later. In Step 2, we show the time decay estimates for the basic energy. 

Step 1. Multiplying (|3.3.1I) by r 3 — x 3 and integrating the resulting equation with respect to the 
spatial variable, we have, with the help of the integration by parts and boundary condition (13.3.2Jh 
that 



Notice that 

55 (r 3 — x 3 ) x — 4Ai (r 3 — x 3 ) = 355 ( r 2 r x — x 2 ) — 4Ai ^ 


v x vr a 


(r 3 - x 3 ) 


=4Aix 2 
=4Aix 2 
which implies 


v v 7 


-h 

r r r r 


xv r xvr 7 


, ( r \ xr x 

In ( -J -I- 

V xr x J r 


+ 3A 2 a: 


+ 3A 2 x 2 


r rv v T v 

~^v x + 2 ~^r x -- - 2 - 

x z r T r 


-rrr x - In —^r x - 1 


J t 


®( r3 - a;3 ),- 4A i il) x ^- xZ ) 


dx = I x 2 
dt 


! 


4Ai<hi 


r 

xr x 


+ 3A 2 $ 2 


dx , 


where 


4>i(z) := Inz + z 1 —1 and $ 2 ( 2 ) '■= z — In z — 1 . 
Notice also that 


(3.3.14) 


x r“ 


Then, we set 


rjo := x fj 0 + x^pv 


r x 


x r* 


and i)o := 4Ai$i - + 3A 2 4 > 2 - 5 r x , 


xr 7 


and obtain 


J t J m(x,t)dx + J fP | ^4 (r 3 - x 3 ) ^ - (jL-'j (r 3 - dx 

= fx 2 pv 2 (1 + 2 ^jdx. 

Noting that the quantity {•} on the left-hand side of (|3.3.16p can be rewritten as 

^y- 3 f^V - (-Yr x + 4(-Yr x -7(-Y + 4 X 

r z r x J \r z r x 


(3.3.15) 


(3.3.16) 
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we can then show, using a similar way as to the derivation of (13.3.71) . that 




> 


3P7-4) f 2g7 


.T p 1 


2 (-1) +(r c -l) 5 


cLc, 


when (13.2.31) holds for a small eo- It follows from (13.3.16[) and (13.3.171) that 


dt 


j r/o(x,t)dx + 


[ x2 fC 

2 (— — l^ 2 + (r x — l) 2 

J 

\x J 


s 


dx<C v 2 dx. 


(3.3.17) 


(3.3.18) 


This implies, with the aid of (13. 3. lip and (I3.2.5p . that 


3(3 7 - 4) 


rj 0 (x,t)dx + —— J 


x p' 




dxds 


<C J(v + Vo)(x,0)dx. 

It remains to analyze r)o- By the Taylor expansion and (13.2.31) . one may get that 


(3.3.19) 


fto>\ 


4Ai ( ^r-l) +3A 2 dr x - 1 


x “ 


3 

- 4 a 


2 1—-l) +(ir x - 1 


xr x 


x- 


where a = min{2Ai/3, A 2 }; and 
Vo < 2 


4Ai I —-1) + 3A 2 (^r x - 1 


xr x 


x‘ 


< C 


G - 1 ) 2 + (r *- 1)S 


(3.3.20) 


provided that eo in (|3.2.31) is suitably small. Notice that 

r \ 2 / r 3 \ 2 / r r 2 A 2 

X _1 ) “V^ 3 ” 1 ) ~ \ ~XTx "x 2 1 ~ C 


\ 2 / 2 

— ~l)+ ( r -r x -l 
xr x J \x z 


< Co 1 % 


and also 


(r x ~ l) 2 < Cfjo + C (^- - lj < Co %. 
We then achieve, with the help of (13.3.191) and (13.3.111) . that 


* I x2 1(1 - + (ri - 1] 


dx + 
+ [r x ~ 


x 2 p 7 


r \ 2 

- - 1 

x 


<37 i] LS 

( r x — l) 2 dxds < C J(rj + rfo)(x,0)dx. 


(3.3.21) 


This, together with ()3.3.8[) and (|3.3.20p . implies (|3.3.121) . 

Step 2. We are ready to show the time decay of the basic energy. Let rj be given by ([3.3.61) . It 
follows from ([3.3. 101) that 

(1 + t) j rj(x,t)dx + 3o j (1 + s) J ^—v x + 2 r x v 2 ^j dxds < J rj(x, 0)dx + J J rj(x, s)dxds. 
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In view of (13.3.81) . (|3.3.3I) and (13. 3. 121) . one has that 


J J r)(x, s)dxds <C J J v 2 dxds + C J J x 2 fP ^— — 1^ + (r x — 1) 


<C 

<C 


dxds 


’ J (vo + v){x,0 )dx. 

So, it holds that 

(1 + t) J r](x,t)dx + 3a J (1 + s) J ^— v 2 + 2r x v 2 ^j dxds < C j (770 4- T])(x, 0)dx. (3.3.22) 
This, together with ()3.3.8[) and (|3.3.20D . implies 

(1 + t) xp^v(-,t) + (37 — 4)(1 + t) xp? (J- - l,r x - l) (-,t) 


+ cr 


fo 


+ s) \\(v,xv x ) (•,s)|| 2 ds < C 


- 1 , ro x - 1 


+ 


xpiv(-,0) 


Since xp 7 1 is equivalent to the distance function, dist(x,dl), it then follows from the Sobolev 
embedding (13.2.81) . ((2.1.311 and the Holder inequality that 


J (r — x) 2 (x , t)dx < J x 2 p 2<a 1 ' ) ((r — x) 2 + (r x — l) 2 ) (x, t)da 


< 


x 2 ((r — x) 2 + (r x — l) 2 ) (x, t)dx 


2-7 

7 


27-2 

<C{l+t) T 


X 


r 0 


- 1 ,r 0x - 1 


+ 


xp 2u(-, 0) 


c 2 fP ((r — x) 2 + (r x — l) 2 ) dx 
2 


27-2 

7 


(3.3.23) 


This finishes the proof of ()3.3.13D . 

□ 

With the estimates obtained so far, we are able to derive a pointwise bound for | r/x — 1| and 
\t x — 1| away from the origin, by realizing that equation (13.3.11) can be integrated with respect to 
both x and t. It should be noted that the monotonicity of the Lane-Emden density which decreases 
in the radial direction outward plays an important role for this estimate. 

Lemma 3.7 Let I 2 = [1/2,1]. For a suitably small constant eo in (|3.2.31) . it holds that for (x,t) G 

h x [0,21, 


x l r(x,t) — 1 + | r x (x,t) — 1| 


<C 


( xp*v{-,ti) + X - 1,7-Ox - l) + II^Ox - 1 |Il°°(/ 2 )) • 


(3.3.24) 


Proof. The proof consists two steps. 

Step 1 (bound for r/x — 1). Notice that 


x[r — x = 


/ X 

\y( r (y, t) — y) 2 ] y dy 


< ||r — x||" + 2 ||r — x[| ||a;(r x — 1)|| . 


This, together with ()3.3.121) . yields that for x € I 2 , 


- 1 


< 8x(r — x) 2 <C ( xpzv(-,fy 


+ 


1 j 1"0x 1 


(3.3.25) 
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Step 2 (bound for r x — 1). Integrating equation (13.3.111 over [x, 1] and using the boundary 
condition (I3.3.2p . one gets 




2 - \ 7 

ar p x 


r 2 r. 


V 


p(-) vtdy-[——\ - -T (p 7 )„ dy = -53 + 4Ai 



rJy 


(3.3.26) 


where 53, defined by (|2. 1.121) . can be rewritten as 


53 = n (In r x ) t - -Ai - 2A 2 (lnr) t . 


So, (13.3.261) is equivalent to 
p( ln r x ) t = f 


x 2 p x 7 
r 2 r T 


1 -V 2 


l -y 2 2 


p—^vdy - 2 / p-rrw dy + 


- 1 P 4 


(p 1 ) v dy 


+ (j A i -2A 2 ^ (lnr) t + 4Ai J (In r) yt dy. 
Integrate it with respect to the temporal variable to obtain 


yin ( — 

\rox 



t / x 2 p \ 7 , 

—o— ) ds + £, 
o \r z r x J 


where 


£ = - 


L 


1 rl 

Pfa vd y 


t ft ri „,2 

— 2 
o 


n pK^v 2 dyds + 


t r V 


+ I — 2A 2 ) In 



0 J x r 
r 


r 0 . 


+ 4Ai In 


o li tr 
r(l,t) r 0 (x) 
r 0 (l) r(x,t) 


(.F) v dyds 


which implies that 


r x = r 0x exp 


—211 exp < — i 

y J [y 


rt / x 2 p \ 7 

where 21 = / ( — I ds. 

r- r. 


(3.3.27) 


On the other hand, direct calculations show, by virtue of (13.3.271) . that 


2h = 


x 2 p x 7 
r 2 r r 


x 2 p x 7 
r 2 r 0x 


P 


exp < ——21 > exp <-£ 


p 


so that 


exp { —21 }> = 1 + 


7 (^_P_\ 


0 p, V r z r 0x ) 


exp \ ——£ [> dr. 


It then follows from (13.3.271) that 


r x =r 0x 


i+ 

Jo y \ r r ox J 


exp \ ——£i )■ exp <! — 


y 


f*T r 1 4 

7 r f y 
y Jo 



(p" / ) v dyds [• dr 


1/7 


exp | —£1 | exp | — [ f 

Ip J Ip7o A 


* 'V 


(y 7 ) dyds ^ , 


(3.3.28) 


where 


£1 =£- 


/ J fl(^%dyds. 

do dx r 
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In view of (13.3.31) and (13.3.251) . one can get that for x > 1/2, 

1/2 


|£i| <c( X 


v 2 pv 2 dx J pdx^j + C x 2 pu^{ro(x))dx J pdx 


1/2 



+ C j j v dyds + C 


- 1 


x 


L°°(J 2 x[0,T]) 


< Cl 


where 


c = 


xpzv(-, 0) 


+ 


/ r ° i i 

X\ - 1 , Tq x — 1 


It therefore follows from (13.3.281) and (13. 3. 251) that 


r x <rox 


1 + 


(1 + Ct) f 
J o 


t T' -7 

—p 1 exp 

o T 

t r i 4 


7 r f 1 y 4 

T Jo 


(. P^)ydyds'j dr 


1/7 


x(l + C C )exp|i/ / 


(1 + Ct) 

U + Ce) [ 
Jo 


exp 
't 


7 f 1 y 4 


(. P y ) v dyds 


+ 1 


7-7 

—p ' exp 
0 


t ri y 4 


(, P 1 ) v dyds \ dr 


1/7 


where e = e + ||ro x — l|li,oo(/ 2 ) • Observe that ( fC) x < 0. So, one can derive from (I3.3.25|) that 


r x <ro x (1 + Ct) 

(1 + Cc) [ 

Jo 


exp 

't 


{U Sc i - c ‘W vis ) 


+ 


7-7 

—p ' exp 
0 h 


7, 


t r 1 ] V7 

(1 - Ce) ( p y ) y dyds dr 


exp -(1 — Ci)p 1 t 

/i 


(1 + Ct) 

+ (1 + C'c) [ — fP exp/ — — (1 — Ct)(P{t — r)l dr 

Jo t { y J 


< r 0x (1 + Ct) 


exp <{ ““(I — Ct)fPt }> + 


1 + Ct 


IJtqx (1 4- Ct) 


| exp | 


--(1 -Ct)^t 

fl 


1 -Ct 
1 + Ct 


1/7 


exp — — (1 — Ct)fP(t — t) 


t 

T —0 


1 1/7 


1 - 


1 -Ct 


+ 


1 + C'e'l 


1/7 


1 - 


Ct j 


< r 0x (1 + Ct) 


Similarly, r x > r$ x (1 — Ct ). These two estimates, together with (I3.3.25p . imply (13.3.241) . 


□ 


The following lemma gives the decay estimates for the weighted norms of both the time and 
spatial derivatives of v. 

Lemma 3.8 Let (13.2.31) and (13.2.41) be true. Then it holds that, for 0 < t < T, 

(1 + t) [xp^v u v,xv^j (-,t) +J (1 + s) \\(vt,xv tx ,) (-,s)|| 2 ds 


< C 


xp2v t ,v,xv x ) (-,0) 


+ 


/ r ° 1 1 

X\ - l,r 0x -l 

x 


(3.3.29) 
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Proof. Multiplying equation (|3.3.1I) by r 2 and differentiating the resulting equation with respect to 
t, we obtain 


px 2 v tt - 7 r 2 


Vf 2 -+^ 

r 2 r. 


+ 2 rv 


2 - \ 7' 

x A p x 

r 2 r v . 


.X 


+ 2 z* v (p r )i 


=r 


Kxt + 4Ai - 


r/xti 


+ 2ru 


®x + 4Ai - 


r / sJ 


Set 


r/i(x,t) :=-x 2 pvf + — — 


2 - \ 7 r 

xp x 


r 2 r 7 


4 r - 3—i-r x v + 2^-uu 


"7 7*^ 

(2q - l)r x u 2 + 2(7 - l)ruu x + u 2 

— r. T 


x 


P 7 


Following the estimates for 7 defined in (13.3.61) . we can show that, for 7 e (4/3,2], 

'v\2 


/ ^ ^ 2- 2 , ^7 4 2-7 

*7i (M) > 2 X P v t + —4 — x p 


x) + V * 


(3.3.30) 


(3.3.31) 


(3.3.32) 


Vi(x,t) < ^x 2 pv 2 + cx 2 p 1 


v\ 


+ < 


x) 


(3.3.33) 


provided that (13.2.31) holds with eo being suitably small. Multiplying (13.3.301) by vt and integrating 
the product with respect to the spatial variable, we have, using the integration by parts and 
boundary condition (|3.3.2I) . that 

( r2 vt) x - 4Ai r 2 v t 


dt 


J rji(x,t)dx + 


v 

r/ xti 


dx = 3i + 3 2 , 


(3.3.34) 


where 


3i := — 


*B(2 rvv t )x ~ 8A1 


r / x 


rvvt 


dx, 


3 2 -VI; - 1)1 




2 - \ 7 
X Z p' 

r 2 tv 


2 - \ 7 2 "l 

x z p \ r z 


r* r« 


' X J t 


v 2 dx + 2(7 — 1) J 

v 2 x dx- f p< 


2 - \ 7 
X Z p X 


r* Try 


VV T dx 


X 


4-7 - 3— T r x V + 2 — vv. 


dx. 


The second term on the left-hand side of (|3.3.34D can be estimated as follows. Notice that 


= fVt 


3 r x \ r / x 


+ A; 


(■ r2y t ) s 


T x r 


+ % 


where 


* : =5 a ‘ 


UN 2 

r> 


Thus, 


( r2 vt) x - 4Air 2 ^ 


— A 2 


v 

r / xti 


2 - + — 


(3.3.35) 


dx 


vt 
r / x 


r 

r x 


(■ r2 vt) x ~ 3r2v t 


+ A ; 


|(r 2 u t ) ; 


r x r 


dx — 33 


(3.3.36) 


>3(7 


v tx + 2 r x v t 


dx — J 3 , 
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where a = min{2Ai/3, A 2 } and 


d_ 

dt 


vtx H - 2 v x Vt I dx < Ji + J 2 + ^3- 


A3 := - J ® (r' 2 ut) a , dx - 4Ai J r 2 v t dx. 

So, (13. 3. 341) implies that 

rn{x,t)dx + 3a 

For 3\ and A 3 , it follows from (I3.2.5|) . (13.2.41) and the Cauchy inequality that 
Ai + A 3 < a j (^—v 2 x + 2 r x v^j dx + Ccr~ 1 ef J (x 2 u 2 + v 2 ) dx. 
Similarly, A2 can be bounded by 

A 2 < Ce 1 J [x 2 v 2 x + v 2 ) dx. 

So, we arrive at the following estimate 

J r)i(x, t)dx + 2a j v 2 x + 2r x v 2 ^ dx < C J [x 2 vl + v 2 )dx, 
provided that (|3.2.4D holds for ei < 1. This, together with (13. 3. Ill) , implies that 

j rji (x, t)dx + a J J [x 2 v 2 x + u 2 ) dxds < J r/i(x,0)dx + C J rj(x,0)dx 


and 


(1 +t) J r/i(x,t)dx + cr J (1 + s) J (x 2 v 2 x + v 2 ) dxds 

r)i(x,0)dx + J J 771 (x,s )dxds + C J (1 + s) J (x 2 v 2 + v 2 ) dxds 
7]i(x, 0)dx + C J J v 2 dxds + C J (1 + s) J (x 2 v 2 + v 2 ) dxds. 
Here (13. 3.33 1) has been used. This, together with (13.3.131) . 


< 


< 


implies 


(1 + t) xp 2 Vt {-,t) + pa (v,xv x ) (-,t) j+J (1 + s) ||(w a ,xt; sa: ,) (-,s)|| 2 ds 

/ 1 2 7 2 /Tc\ \ 2\ 

< c ( xp2(v t ,v)(-,0) + p 2 (v,xv x ) (-,0) + 1,7-Ox -lj )• 


Observe that 


(1 + t)v 2 (x, t) < u 2 (x, 0) + f (1 + s) [2 u 2 (x, s) + v 2 (x, s)] ds. 

Jo 


1 0 

Integrate the above inequality with respect to the spatial variable to give 

(1 + 1) / v 2 (x,t)dx < J v 2 (x,0)dx + J (1 + s) J [2 v 2 (x, s) + v 2 (x, s)~\ dxds. 


(3.3.37) 


(3.3.38) 


(3.3.39) 
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Similarly, it holds that 

(1 + t) J x 2 v 2 (x,t)dx < J x 2 v 2 (x,0)dx + J (1 + s) J [2x 2 v 2 (x, s) + x 2 v 2 (x, s)] dxds. 

This, together with A3.3.131) and (1-1. -1. ■101) . implies that 

(l + i)f xp?v t (-,t) + \\{v,xv x ) (-, t)|| 2> ) + f (1 + s) ||(u s ,am sa; ,) (-, s )|| 2 (is 

V ' Jo (3.3.40) 


< C 


xp 2 v t (-, 0 ) + \\{v,xv x ) (-, 0)|| 2 + 


, r ° i i 

x \ 1 , rox 1 

x 


This finishes the proof of ({3.3.29(1 . 


□ 


Next, we derive further time decay estimates based on Lemmas 13.5113.61 and 13.81 bv using two 
multipliers 

rx nx 

/ P _/ 3 (y)(r 3 - y 3 ) y dy and / p~ p (y)(r 2 v) y dy for 0 < /3 < -y — 1 . 

Jo Jo 

The key is to deal with the behavior of solutions near both the boundary and geometrical singu¬ 
larity at the origin simultaneously. The improved decay estimates obtained in this lemma give the 
convergence of the evolving boundary r = R(t) to that of the Lane-Emden stationary solution. 


Lemma 3.9 Suppose that (|3.2.3D and (13.2.41) hold. Then for any 6 £ (0, 2(7 — l)/(37)), there 
exists a constant C{6) independent oft such that 

7$ 7—1 2 7~ 1 n g 3(7 — 1) r \ 

p~ ~ (r-x,xr x -x)(-,t) + (1 +t) y || (xr x - x) (-, t)\\~ + (1 + t) r 

x ||(r - x)(-,t)\\ 2 + (1 + t )~i (xp5v t ,v,xv x 'j (-,t) + p% (r - x,xr x -x)(-,t) ^ 


+ 


+ 


L 

L 


07+2 


p 4 (r — x, xr x — x) (•, s) + (1 + s) 


7~1 


-e 


pa (r — x, xr x — x) (-, s) 


ds 


(1 + «)" 


27-1 

7 


-0 


27-1 e 


(v,XV x ,V s ,XV sx ) (-, s) 11 *^ + (1 + s) 2 r 2 P 4 2 (v,XV x )(-,s) 


7 6 _ 7 —1 


1 ]) 


< ( 7 ( 0 ) ^ (v,xv x ,xp 2 v t ^ (-,0) + Wox ~ l|li°o^ , t £ [0, T], 

Moreover, we have for any a £ (0,1) and 6 £ (0, 2(7 — l)/( 37 )), 

(1 + t) 2(7-D/7-»|| (r _ x)(.,t)||J,„ ([a l|) + (1 + t)(27-l)/7-»||„(.,f)||| 00(h 

<C(a,0)f pxv x ,xp 1 ^ 2 v t 'j (■, 0) + IIt'Ot. - 111!™) , te[ 0 ,T]. 

Proof. For any given 6 £ (0, 2(7 — l)/( 37 )), we set 

„ 7 l0 6 p 1 /. £ \ 

/3:=7-l- T . ‘“j, 1/: =2( 1 + 7~V' 


ds 
(3.3.41) 


(3.3.42) 


(3.3.43) 


so that 0 < P < 7 — 1 and 0 < l < P/(2j). The proof of this lemma consists of the following four 
steps. 
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Step 1. In this step, we prove that 

/ [(1 + W + 1]*V' - l) + (r x - l) 2 ( x,t)dx 

+ j (1 + S Y J P~P { x2y l + v 2 ) dxds + J fx 2 p^ - l) +(r x -iy 
< c \\r 0x - lllioc + C^j\l + 8y\Ki\da + C Y.[ AA, 


dxds 


(3.3.44) 


where 


L i = - J P^vt ^j P P(r 3 - y 3 ) y dy \ dx, 

p~P ( r 3 _ _ f p-l 3( r 3 _ y ^ ydy 

Jo 

[ p~P(r 3 - y 3 ) y dy - p~P (r 3 - x 3 ) 


L 2 = <? -j 



L 3 —4Ai 



v 
r/ x 


dx, 

dx\ 


(3.3.45) 


and 


K\ =- j p^v t (/ p P (r 2 v) y dy S J dx, 


K 2 = j fP 
K :i =4Ai 



x 
r ~4 

V 
rJ x 


p 13 r 2 v 


[ P P {r 2 v) y dy 
Jo 


dx, 


(3.3.46) 


p l3 (r 2 v) y dy ) — p ^r 2 " 


dx. 


To this end, we multiply (13.3.11) by fj p~^(y)(r 3 — y 3 ) y dy and integrate the resulting equation 
with respect to the spatial variable to obtain, with the aid of the integration by parts and the 
boundary condition (I3.3.2p . that 


+ J fr 

Noticing that 


/{ [£ (r 3 - * 3 )] _ - (^r) 7 (r 3 - * 3 ) x } * 

®{r 3 -x 3 ) x - 4\! (I) {r 3 -x 3 )]dx = Y J L l . 


2—1 


J x 2 p 13 (1 - 1) + (r x - 1 y 


(x, 0 )dx 


<C\\{r 0 -x, xr 0x - x)|| 2 oo f (1-x) A dx < C ^ 0x 

Jo ( 7-1 )-P 

due to (|2.1.3D and ro(0) = 0, one can obtain, following the derivation of (13.3.211) . that 


x 2 pP 


Gr 1 ) 2 + (r *“ 1)i 



(x, t)dx + / x 2 p 1 13 


-1 +{r x -iy 

x 


<C\\r 0x -l\\ 2 Loo +C 


E f a a 

i=l J ° 


dxds 

(3.3.47) 
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Next, multiplying equation (13.3.11) by f^p ^{y){r 2 v) v dy, integrating the product with respect to 
spatial variable, and using the integration by parts and the boundary condition (|3.3.2j) . one obtains, 


p<-P 


v - 


„2 \ 7 


r~r n 


( r 2 v) x dx + I p 13 53 {r 2 v) x dx — 4Air 2 u dx = ^ Kj, 


i= 1 


Following the derivation of (13.3. 10p . one can then obtain 

^ J r] 2 (x,t)dx + 3n J p~ f3 


fr" 

—v 2 x + 2 r x v 2 

J 

l r x 


dx < 22 Ki, 


1=1 


where 


r] 2 (x,t) := p & I rj(x,t) - ^x 2 pv 


1 


2 - ,2 


I 2 ^ f 


G _1 ) +(r * _ir 


Here and thereafter, / ~ g means that C~ l g < / < Cg with a generic positive constant C. 
Multiplying the equation above by (1 +t) v and integrating the product with respect to the temporal 
variable lead to 

(l + t) 11 J x 2 p y ~ 13 - 1 ^ + (r x — l ) 2 (x,t)dx + J (1 + s) u J p~P (x 2 v 2 + u 2 ) dxds 

< C \\r 0x - l\\ 2 Loo + C 22 j (1 + s) v \Ki\ds + cj J x 2 p 1 ~ 13 - l) + (r x - l ) 2 dxds, 

(3.3.48) 

due to the fact @/7 < 1. So, estimate (|3.3.44p follows by a suitable combination of (I3.3.48D and 
(13.3.471) . 

Step 2. In this step, we show that 


(1 + t) 2u pi + (1 + t) K + p-P 


x 


P-l 2 + (r,^f 


(x, t)dx 


x 2 fp 




+ (1 + t) 2u J x 2 pv 2 {x,t)dx + J J (1 + s) k + /9 13 

+ j j (1 + s) 2v + (1 + s) v p-P ( x 2 v 2 + u 2 ) dxds 

<c( xp?v(-, 0) + ||r 0x - lUioc^) + C 22 [ + s) u \Ki\ds + C22 ( \ L M S , 

V / i=i Jo i=i Jo 


dxds 


(3.3.49) 


where Li and fQ ( i = 1,2, 3) are given by ()3.3.45p and (|3.3.46p . respectively. 

To prove (|3.3.49p . one can integrate the product of (1 + t) 2u and (13.3.101) with respect to the 
temporal variable to get 

(1 + t) 2v J ^x 2 pv 2 + x 2 f (J- — l) +(r x — l ) 2 | (x, t)dx 
+ J (l + s) 2u J {x 2 v 2 + v 2 ) dxds < C xp2v(-, 0 ) + C \\r 0x - l\\ 2 Loo 


/ 0 

+ C J (1 + s) [ v 2 dxds + C 


J v 2 dxds + C J (1 + s) 


2-7 
X p' 


G - 1 ) 2 + (r *- 1): 


dxds, 
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since (3 /7 < 1. Integrate the product of (1 + t) K and (13.3.181) with respect to the temporal variable 
to give 


(1 + t) K J x 2 — lj +(r x — l) 2 ( x,t)dx 
+ j\l + sY J x 2 fP g -iy + (r x -iy 


dxds < C || r 0x - 1\\ 2 lo 


+ C J (1 + s) J v 2 dxds + C J (1 + s) L 1 J x 2 (J- — l) + {r x -iy 


(x, t)dxds. 


The last term on the right-hand side of the inequality above is estimated as follows. It follows from 
the Holder inequality and the Young inequality that 


/' 


- 1 + (r x - iy 

x 


dx<[ x p 




X 


--1 Y + (r x -iy 


dx 


x I / x 2 p & 


G _1 ) +(r * _ir 


dx 


1~P 

7 


and 


j\l + s fh-^' Jx 2 (I - l) 2 + (r„ - 1) 

(F 1 )’ + (r ^ 1) 


dxds 



x 2 p 7 P 


< C 


< C1 —£ sup [ x 2 p~P 


dxds 


f (1 + s)^ 7 1 sup f x 2 p P 

[(^-l) 2 + (r,-l)4 

Jo sG[0,t] J 



(x,s)dx 


+ C 


L 1 sG[0,t]. 

x 2 fP~P 


Gr 1 ) + (r * _ir 



/ 0 


(x, s)dx 


dxds. 


y- 1 ) +(r *- i): 

In a similar way as to deriving ()3.3.44|) . we then have, noting ([3.3.131) . that 
(1 + t) 2u J l^x 2 pv 2 + x 2 p 1 - lj +(r x - l) 2 | (x,t)dx 


+ (1 + t) K J x 


— lj + (r x — l) 2 (x,t)dx + J (! + .)*/ (x 2 v 2 + v 2 ) dxds 


< C 


/ 

x 2 fp 

- l) 2 + (r* - l) 2 j 

J 


Vx / 

0 ) 

+ C II r Qx - 11| rex, + sup / 



sG[0,t] J 


dxds 


(3.3.50) 


(x,s)dx 


+ C 


a 


x2 p r p +(^-i ) 5 


dxds. 


Make a summation of A; x (13. 3. 441) and (13.3.501) with suitable large k to give (13.3.491) . 
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Step 3. We claim that 

[[a +t)^- e f+{i+t)^- e +p^- l -^ 6 )]x 2 g-i ) 2 + (r,-ir 


(x, t)dx 


+ j J (1 + s) y 6 + (1 + s) 2 t 2 p ( 7 1 270 ) (x 2 u 2 + u 2 ) dxds 


+ 



x 2 p 7 


G _1 ) +(rs_ir 


’(l + s)V- e + p-( 7 - 1 -l 70 ) 

+ (1 + i)~ 0 J x 2 pv 2 (x, t)dx < CQ( 0 ), 

where and in the following 

Q( 0 ) := (v,xv x ,xp^v t ^j (■, 0) + ||r 0a; - l|| 2 oo • 


dxds 


(3.3.51) 


To prove this claim, it remains to estimate LQ and L* in ([3.3.49(1 . First, it follows from (12.1.31) 
that for any given constants 5 £ ( 0 , 1 ] and /3 £ ( 0,7 — 1 ), 


L{I« r ' 3 dy ) dx - c Li 1 -^ t' 


1 — (1 — x ) 1 T - 1 


dx 




(3.3.52) 


Let oj £ (0,1/2) be a small constant to be determined at the end of this step. It follows from the 
Cauchy inequality, (13.3.291) . the Holder inequality and (13. 3. 5211 that 


J (1 + sY\Ki\ds < u J (1 + S Y J P ^(M + \yv y \)dy^ dxds 

+ Coj~ 1 f (1 + s) f v 2 dxds < Cuj~ 1 Q(0) + Coj f (l + s) 1 ' f p~^(\v\ 2 + \yv y \ 2 )dyds 

Jo J Jo Jo 


(3.3.53) 


since 



p P y(\v\ + \yvy\)dy^ dx < J p ,3 (\v\ 2 + \yv y \ 2 )dy J (^J P P y 2 dy)dx 


<C 


[ P P (\v\ 2 + \yv v \ 2 )dy. 

Jo 


Similarly, one can obtain 

f \L\\ds < Cuj- 1 f f v 2 dxds + ui f / p 2 ( / p~^y 2 

Jo Jo J 


< Cuj-'Qi 0) + Cuj x 2 p 7_/3 



/ 0 j vo 
2 


- - 1 

y 



+ \r y — 1 | dy ) dxds 


Gr 1 ) +( r *~ 1 Y 


(3.3.54) 


dxds, 


since 


J P 2 (^J P 7 dx < C J p 2 p 7 d+Pi l ^dx = C J p 2 7 p /3+ ^ 7 [S> dx < C. 
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X 2 can be rewritten as 



p /5 ) r 2 vdydx 
' y 

f p~P (r 2 v) y dy 
Jo 


p Pr 2 v 


dx =: K 2 \ + K 2 2 ■ 


Note that 

1^211 = 


P f 1 -^ 

~ / I 

7 Jo 


P 7 



p P ( 7 1 Jy(f>r 2 vdydx 

\ 1/2 


<C £ “ p^x( K \r x -l\ + \ r --l[)x- 2 (^j* y Q dy^ ' dx(^j\ 2 dy^j ' 


<c(^Jx 2 p' y (J- - i) + (r x - l) 2 ^ v 2 dy^j 

<C{1 + t)~ v J x 2 ff — 1^ + (r x — l) 2 dx + Coj( 1 + t) v J v 2 dx 


1/2 


and 


l^22| < P 1 -4 


' 1 —UJ 



p Pr 2 v — p P(r 2 v) y dy 

Jo 


dx 


-I 


< C J xp 1 P (jr x - 1| + 

r 1 


+ C J xp 7 (Jr x - 1| + 

(1 + ()-'£ 
0 + t)-'I 


r 

- - 1 

x 


- 1 



\v\dx 


P 13 (M + | yv v \)dy ) dx 


~ —UJ 

1 


x 2 (P~P 




x 2 (P~P 


\r x ~ 1| 2 + 


\ r x — 1| 2 + 


-1 


- -1 

X 


v 2 dx 


dx + (1 + t) u J 

dx + (1 + t) v f ( v 2 + x 2 v 2 )d: 

Jo 


due to 


fi P * +l3 (/ P + \ yVy ^ dy ) dx 

<( [ P^ +l3 [ p~ 2 P dydx\ [ (M 2 + \yv y \ 2 )dy <C [ (|u| 2 + \yv y \ 2 )dy. 

\Jl-uj JO J Jo Jo 

Then, one gets, using (|3.3.121) . that 

J (1 + sY\K 2 \ds < Cuj~ 1 Q( 0) + C (u + + W 2 fRij 


Uo 


(\r x -l\ 2 + 

- - 1 

2\ 

dxds 

V 

X 

J J 


(3.3.55) 


Similarly, L 2 can be rewritten as 



P y3 ) ( r 3 - y 3 )dydx 
J y 

p~P{r 3 - x 3 ) - [ p~P{r 3 - y 3 ) y dy 
Jo 


dx L 21 + L 2 2 - 
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Note that 


1^211 - 


p_ r 1 -" 

7 do 


p 7 



<c 


I fpx (|r x — 1| + |^ — l|) 


p & ( 7 1 ^y(/)(r 3 — y 3 )dydx 

p-P~('V~ 1 )y 2 \ r — y\dy ] dx 


—2 


<cv) J 


2 -7 


--l ) 2 + (r a; -l ) 5 

x / 


dx 


and 


|-^ 22 1 <c f xp 1 13 (\r x - 1| +-1 ) |r - x|dx 

dl-w v x J 


+ C xp 7 (|r x - 1| + 


- 1 





<w J x 2 pP~ p 


p‘ 


p V ( |r„ - 1| + 

i 



- 1 


dy ) dx 


+ C( x z p 


r 2 7f1~P 


pl+P 


' l—o; 


f-l) + (r x — l) 2 dx + Cu 1 [ p 1 /3 (r —x) 2 dx 

Kx 1 J Ji- U 

/ r \ 2 2 1 \ 1/2 
(- - !j + (r x ~ 1) dx) 

^ dy'j dx J 


p V ( l r ?/ “ !| + 


2 \ 1/2 

V 


<C f (up 1 13 + u V 7 ) z 2 _ 1 ) + ( r z “ 1 )" 


dx, 


where we have used the following simple estimates due to (13.2.10p and (12.1.31) : 
f 1 (F~P{r - xfdx <C f 1 (P~P+ 2 (^-i) [(r - x ) 2 + (r* - l) 2 ] 

J 1 - 0 ) J1/2 

’ [ f [(r -x) 2 + x 2 (r x - l) 2 ] 

J 1/2 


<c 


and 


r 

- - 1 
y 



dy ) dx 


< C 


< 


P /P+e f p~ 2J3 ~ 1 {y)dy [‘ iPy 2 ( |r„ - 1| 2 + 
di-w do do y 

C f p~^dx f x 2 p 1 (\r x — 1| 2 + 
dl-w do V 

< c J x 2 p 7 f|r x - 1| 2 + 


- 1 


dydx 


- - 1 

x 


Thus, it follows from these and (13.3.121) that 


[ \L 2 \ds < C{u)Q{ 0) + Cu [ [ x 2 pi-P 

Jo Jo J 


I r x - 1| 2 + 


- 1 


dxds. 


(3.3.56) 
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Rewrite K 3 as 


K 3 = - 


4X1 L G)J„ {^) y (r 2 v)dydx 

/ 


— 4Ai 

1 -u Vr/ z 

K 3 \ can be bounded by 


p ^r 2 v 


P P {r 2 v) y dy 


dx =: K m + K 32 . 


\K 3 i\ = 


<C 

<c 

and K 3 2 can be bounded by 


4—Ai / [ P 13 ^ 1 ' > y<f>r 2 vdydx 

7 Jo x Jo 

J (\xv x \ + \v\)x~ 2 (J y 3 \v\dy\ dx 

/ l —UJ / rl \ 1/2 n 

{\xv x \ + H) yj v 2 dyj dx < C J {y 2 + x 2 v 2 )dx, 


\K 32 1 <c [ p P (\xv x \ + \v\)\v\dx+ [ {\xv x \ + \v\)( f p 13 (|v| + \yv y \) dy] dx 
Jl-u Ji-uj \J 0 J 

<u> f p~^(\xv x \ 2 + \v\ 2 )dx + w _1 f p~ l3 \v\ 2 dx 

J 1—U) J 1— UJ 

+ uj~ 3 [ (\xv x \ 2 + \v\ 2 )dx + UJ f p~ 13 (\v\ 2 + \yv y \ 2 )dy, 

Ji-lj Jo 

due to ()3.3.52p . Since /3 < 7 — 1, the Hardy inequality (|3.2.10l) implies that 

f p~P\v 3 dx < C j p~ /3+2 ^' y ~ 1 \v 2 + v 2 )dx < C f (v 2 + v 2 )dx. 

Jl-oj J 1/2 J 

These, together with (13.3.131) . yield 

f\i + s y\K 3 \ds 
Jo 

<Cw _1 f (1 + s) j \x 2 v 2 + v 2 )dxds + Cuj f (1 + s) u f p~ l3 (\xv x \ 2 + \v\ 2 )dxds (3.3.57) 

Jo J Jo J l-uj 

<C(uj)Q(0) + Cuj f (1 + s)" f p~ ,3 {\xv x \ 2 + \v\ 2 )dxds. 

Jo Jl-u 

Similarly, L 3 is rewritten as 

"1 —UJ 


L3 — — 


“ Al i Oj„ (p-^y-y^dyi. 


— 4Ai 


/ l-w Kr/ X 

Clearly, L 3 \ and L 32 can be bounded by 


p 13 (r 3 -x 3 ) - ( j p /3 ( r 3 - y 3 ) v dy 


dx —: L 31 + L 32 . 


\L 3 i I < C J (y 2 + x 2 v 2 )dx + C J p 7 (r — x) 2 dx 
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and 


| £ 32 1 < C / p 13 {\xv x \ + \v\) \r - x\dx + C f (\xv x \ + \v\) f f p Py\r y - l\dy\ dx. 
J 1-0) Jl-O) \ 7 o / 


The second term in L 32 is bounded by 


' 1—CJ 


( \xv x I + M) ( / p ^vVy - l \ d v) dx 


0 


<Cuj 2 


(1 + t) 2u f (\xv x \ 2 + \v\ 2 )dx 
Jl-U 

+ Cuj~^(l + t)~ l ~^ +b ( p~Py 2 (r y -l) 2 dy 
Jo 


1 px 


d dydx 


l-u Jo 


<Cuji 


(1 + t) 2u f (\xv x \ 2 + \v\ 2 )dx + (1 + t) 1 r +t f P f3 y 2 (r y - l) 2 dy 

J 1 —Ld JO 


due to (13.3.521) . The first term in L 32 can be bounded as follows. Since /3 < 7 — 1, it follows from 
(12.1.31) and the Hardy inequality (13.2.101) that for 7 > 4/3, 


P P i\xv x \ + |u|) |r — x\dx 


f 

J1-1 

h f rl r l 

<Cu 2 ^-P (1 + t) u p~^ (v 2 + v 2 ) dx + (1 + t)~ u / p~ d ~ h \r — x\ 2 dx 

J 1 —UJ J l—Ld 

(l + t) u f p~ d {v 2 + v 2 ) dx + f p y ~' 3 [(r - x) 2 + x 2 (r x - l) 2 ] 
Jl-UJ J 1/2 


<C'w 2 ^- 1 > 


dx 


+ (l+f) h + 2 —y f p P [( r _ x ) 2 _)_ x 2 {r x — l) 2 ] dx 
J 1/2 


'1/2 

where h = min{/3/8, (7 — 1 — /3)/4} (it should be noted that /3 + h < 7 — 1), and we have used 
the estimate 

p- 0 ~ h \r - x\ 2 dx < f 1 p 2 ^~ l) ~ 0 ~ h [|r - x\ 2 + (r x - l) 2 ] dx 
. 7 1/2 

h+2 —7 


' 1—UJ 

( rl 


< 


1/2 


p & [(?— x) 2 + x 2 {jr x — l) 2 ] dx 


1 1/2 


p~i P [( r _ x ) 2 + x 2 {r x — l) 2 ] dx 


2 ( 7 -l )-h 
7 


Consequently, taking into account of 13.3.31) . (|3.3.12|) and 13.3.521) . one gets that 


\L 3 \ds <C(uj)Q(0) + C 


r rt 


UJ 2 + (U 2 ^- 1 ) 



(1 + s) 2v + (1 + s) v p- 


x{\xv x \ 2 + \v\ 2 )dxds + sup / p 13 [(r — x) 2 + x 2 {r x — l) 2 ] dx 

[o,t] 7 

+ J J p y ~ 13 [(r — x) 2 + x 2 (r x — l) 2 ] dx 


(3.3.58) 
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We finally derive from (|3.3.49l) and (|3.3.53I) - (I3.3.58I) . by choosing u £ (0,1/2) suitably small, that 


(i+t) 2 v p / + (i + t) K +r fi 


X 


X 


--1 +(r x -iy 


(x,t)dx 


+ (l + t) 2u J x 2 pv 2 (x, t)dx + J J (1 + s) K + p 13 x 2 p y + (r x - 1) 

+ J J (1 + s) 2v + (1 + s) u p~P (x 2 v 2 + v 2 ) dxds < CQ(0). 


dxds 


Due to (13. 3.43 II . this completes the proof of (13.3.511) . 


27-1 

7 

t 


Step 4■ Multiply equation (13.3.371) by (1 + 1) 9 and integrate the product to deduce that 


(1 + f) f 9 J [pv 2 + p 7 (y 2 + x 2 v 2 y\ dx + J (1 + s) 9 J (x 2 v 2 x + v 2 ) dxds < CQ(0 ) 

In a similar way to the derivation of (13. 3.23 p and ()3.3.40p . one can show 

.3(7-!)+g 


and 


f (r — x) 2 (x, t)dx < C( 1 + 1)" 
J (y 2 + x 2 v 2 ) (x, t)dx < C{ 1 + t) 


V( 0 ) 


2 7~1 1 a 

— + Q( 0 ). 


This finishes the proof of (13.3.411) . 

Moreover, it follows from (I3.3.4ip that for x £ [0,1], 

2 (x,t) = (yv 2 (y,t)) y dy < 2 (^J v 2 {y,t)dy^J y 2 v 2 y (y,t)dy 


xv 


+ J v 2 (y,t)dy < C(l + t)~ 


27-1 


r +e Q{ 0). 


Similarly, 

This finish the proof of (13.3.421) . 


x (r(x, t) — x) 2 < C{1 + t) r +9 Q(0 ). 


(3.3.59) 


□ 


3.4 Higher-order estimates 

In this subsection, we derive the higher-order part of the a priori estimates for the strong solution 
(r,v) on the time interval [0, T] defined in Definition 1 2. II under the assumption (13.2.31) and (13.2. 41) . 
To obtain the higher-order estimates, we define 

Q := lnr x + 21n . (3.4.1) 


(3.4.2) 


This transformation between Q and r is one-to-one, and we can solve for r in terms of Q by 
f r x \ V3 

r(x,t) = (3 J y 2 exp {Q (y, f)} dyJ for xGl and t> 0. 


Indeed, we will show in Section 13.4.11 that Q ~ r x — 1, Qt ~ v x , Q x ~ f xx and Qtx ~ v xx . Then 
equation (12.1. llal) can be written in the form of 


x 2 pY x 2 _ 

^ + 7 ( J G*= ~ 


x 


2 \ 7 


r 2 r 7 


x4>p. 


(3.4.3) 


(This is the same as (13.2.71) . we recall it here for the convenience of readers). 
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3.4.1 Preliminaries for higher-order estimates 

The main goal of this subsubsection is to derive some preliminary estimates for the strong solution 
(r, v) on the time interval [0, T] defined in Definition 12.11 and to prove the equivalence of the 
functionals (£(t) and £(t), under the a priori assumptions (13.2.311 and ()3.2.4I) . 

We illustrate how to use Q and its derivatives to control r and v and their derivatives by 
identifying the principal parts of Q, Q t , Q x and Q x t- Note that 


Q = (r x ~ 


1)+2 (i _1 ) +0 ( l 


r x - ir+ 



(3.4.4) 


„ r xx x fr' 

Qx = — + 2- 

r x r V x / x 


fn 4 2 


r 

xJ x i 


+ 


— 1 M XX + 2 (-1 


XJ x 


Qt — -b 2— — ( v x + 2— ) + (- 1 ) v x 2 [ -1 


(3.4.5) 

(3.4.6) 


Qxt = 


V x V 

— + 2 - 

r x r 


— ( Vx + 2 - 


x/X 


+ 


~ 1 1 v xx + 2 (-1 


x 


V 

XS X 


x\- 


~^v x + 2 - - - 


.r/ \x/ x X 

Thus, it follows from (12.2.2l) - (|2.2.4p that for t € [0, T], 

(r x , r/x, v x , v/x, Q, Q t )£L°°(I), (r/x, v/x, v/r) G H l (I), Q x t€L 2 (I), 
(r, v) G H 2 ([0, a]) for a € (0,1), 

(r, v) G # 2 (/), if G L 2 . 

In fact, (13.4.101) can be derived easily by noting that 

and Vxx — Vx Qxt y Vxx 2 ( j 
rt \r / x 




X 

(r\ ‘ 

r XX — fX 

Qx 

- 2- 
r 

\xJ x- 


With the regularity (I3.4.8D - (I3.4.10I) . we have the following Lemmas. 


(3.4.7) 


(3.4.8) 

(3.4.9) 
(3.4.10) 


Lemma 3.10 Suppose that (|3.2.3p holds for a suitably small £q. Then for t G [0, T], 


\\(v x ,v/x)\\ 2 <4\\Q t \\ 2 , 

(3.4.11) 

\\(r x -l,r/x-l)\\ 2 <4\\Q\\ 2 , 

(3.4.12) 

1 \(r xx , (r/x) x )\\ 2 < 4||^|| 2 , 

(3.4.13) 

||(^xx, (v/x) x )\\ 2 < 4 H^tell 2 + c| \(v x ,v/x)\\ 2 Loo \\Q x \\ 2 . 

Here the estimates (|3.4.13p and (|3.4.14p hold if ||<7 r | < oo. 

(3.4.14) 
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Proof. The proof consists of three steps. 

Step 1. In this step, we prove (I3.4.11j) and (13.4.121) . Let s £ (0,1/4) be an arbitrary constant, 
and Xe G [0,1] be a cut-off function satisfying 


Xe = 0 on [0, e] U [1 - 1], Xe = 1 on [2s, 1 - 2s], 

Xe = x/s — 1 on [s,2s], Xe = if ~ x)/s — 1 on [1 — 2s, 1— e]. 

Note that v x = x(y/x) x + v/x and v/x £ H 1 , we use the integration by parts to get 


(3.4.15) 


/ Xc v x v -ix = / X, Q 2 dx + 1 /(0 !2 ) dx = 1 /0 !2 dx - 1 /Q 2 dx. 


Thus, we have 


Xe 


ftr + 2 — 
X 


dx = Xe 


J 


■Z A I V ' 2 

vt + 4 ( - 
x 


I 


dx + 4 / XsVx-dx 


J Xe v‘x + 6 (-^) dx-2jx' £ x(^j dx 

J X , 4 + 60 


-/ X ' 

Note that 


^ + 6 


dx — 2 


dx — 4 


x fv\ 2 , 

— ( — ) dx + 2 


£ VX/ 


r»l—e 


'l-2e 


- f-V dx 


£ Vx/ 


r 2e 


dx. 


v x + 2- = G t - --1 U —2 --1 - 


Thus, 


Xe 


+6 


dx < Xe 


dx + 4 


L 


2e 'V\2 


<2 / Xel/t 2 dx + C6q J Xe 
which implies that 

J Xe 0 + y 0 dx<2j Qlix + 4 J 0 ) dx, 


* + <i) 


dx + 4 


-) 

x/ 


e 
2 

) dx 


dx 


x/ 


(3.4.16) 


provided that (|3.2.3I) holds for a suitably small number eo- Since (v x ,v/x) £ L 2 due to (I3.4.8p . we 
can obtain by letting £ —>■ 0 and using the dominated convergence theorem that 


1 2 11 /u\2 

2*- + y(i 




dx <2 j Q 2 dx + 4 


-Y dx, 
x/ 


which means 


+ 


dx < 4 J Q 2 dx. 

This finishes the proof of ()3.4.11l) . Clearly, (|3.4.12l) follows from similar arguments. 


un2 
x 
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Step 2. In this step, we prove (13.4.131) under the assumption \\G x \\ < oo. Due to Q x £ L 2 (I) 
and (13.4.101) . we have r £ H 2 (I). Then, a similar argument to the proof of (13.4. lip will lead to 
the proof of (13.4.131) . which goes as follows. For the cut-off function Xe defined in (13.4.151) . we use 
r £ H 2 (I) which ensures (r/x) x £ 1]), and the integration by parts to get 


/ 


Xe r xx 


dx = 


fx,( 

!/* 


r 


) (-) dx = f Xe 

x) XX \x) X J 

r 


\lx-\ 


XX £ 


which implies 




dx 



This, together with (13.4.51) and (13. 2. 3D . gives that for small eo in (13. 2. 3D . 


/ 


Xe 



(3.4.17) 


Therefore, by virtue of (I3.4.8D and (|3.4.10D (which implies (r xx , (r/x) x ) £ L 2 ), we obtain (13.4.13D 
with the help of the dominated convergence theorem. 

Step 3. In this step, we prove (|3.4.14D under the assumption ||£7 X || < oo. In view of (I3.4.7D . it 
follows from the Cauchy inequality and (13. 2. 3D that 




X 


2 


cel 






Due to G x £ L 2 {I) and (13.4. 10p . we have v £ H 2 (I). So, we can use the same way as that for the 
derivation of (13.4.17D to obtain 


Xe \ -jV xx + — 


^2 / Gtx.dx + 4 


v 

XJ X 


dx 


-I 


V 

XS X 


dx + c 


v 

v x , - 

X 



(3.4.18) 


where Xe is the cut-off function defined in (13.4. 15D . With the aid of (I3.4.8D and (13.4.10p , we see that 
all the quantities appearing on the right-hand side of (13.4.18p are finite. Since (v xx , (v/x) x ) £ L 2 
due to (|3.4.8p and (|3.4.10D . we can get, by letting e —> 0 and using the dominated convergence 
theorem, that 


/ 



dx < 4 


G xt dx + c 




(3.4.19) 


This, together with (|3.4.13p . gives (I3.4.14p . 

□ 
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Remark 3.11 In view of (13.4.81) . we know that G x t E L 2 {I), t E [0, T], However, this does not 
mean that Q x (t) E L 2 (I), t E [0, T\, unless we assume that Q x (0) E L 2 (I), because 

Gx(t) = Gx( 0) + [ Gxs(x,s)ds, t E [0, T\. 

Jo 

G x { 0) E If 2 {I) is an additional regularity assumption of the initial data other than l£(0) < oo. By 
(13.4.131) . (13.4.51) and (|2.2.13p . we know that the condition Q x ( 0) E L 2 (I) is equivalent to r xx {x, 0) E 
L 2 (I ) if <£( 0) is small. 

Remark 3.12 We can see from the proof of (|3. 4.141) . in particular, (13.4. 191) . that 

\\{v xx , {v/x) x )\\ 2 <4\\Gtx\\ 2 + c\\(v x ,v/x)\\ 2 L oo\\{r xx , (r/x ) x )|| 2 , (3.4.20) 

if \\r xx \\ < oo. To bound \\{y xx ,(v/x) x )\\, we have to control the product of \\(r xx ,{r/x) x )\\ and 
\\{v x ,v / x)\\l<=° , in addition to the bound of ||t/xt|l- Indeed, we prove in Lemma \3. 2f \ that although 
||(r xx , (r/x) x )(-,t)\\ may grow with respect to time, \\(v x ,v/x)(-,t)\\Loo decays faster so that the 
product is bounded. 


Remark 3.13 Similar to (13.4. lip , we can obtain that for any a E (0,1), 

IK^XXI i. v / x ) X )\\ 2 L 2 ([a,l]) - 4 \\Gtx lll / 2 ([a,l]) + C II ( v x, V / x) || Loc([ al ]) \\Gx 11 Z,2 ([ a>1 ]) > 
which implies 

a 2 \\(v xx , («/x)x)||i a([0jl]) < 4\\xGtx\\ 2 + c\\(xv x ,v)\\ 2 LOO \\G X \\ 2 , (3.4.21) 

provided that \\G X \\ < oo. 

Lemma 3.14 Let 5 > 0 be any fixed constant. Suppose that ()3.2.3p holds for a suitable small eo, 
then for t E [0, T], 


P S (r xx , (r/x) x ) 


+ 


xp S ^ 1 ^ 2 {r/x) 3 


< c 


p 5 G* 


P 5 iv xx ,(v/x) x ) 


< c 


p s G xt 


+ c\\(v x ,v/x)\\ z LX , p S Gz 


P 5 (r x ~ l,r/x - 1) 


+ 


/-(7-i)/2 ( r _ x ) 


< c 


p 5 G 


p s (v x ,v/x) 


+ 




< c 


p s Gt 


(3.4.22) 

(3.4.23) 

(3.4.24) 

(3.4.25) 


Here the estimates (13.4.221) and (13.4.231) hold if ||//t? x || < oo 


Proof. The idea of the proof of this lemma is similar to that of Lemma 13.101 Due to p{ 1) = 0, we 
only need to cut off the origin. Let e E (0,1/4) be an arbitrary constant, and rj e E [0,1] be a cut-off 
function satisfying 

% = 0 on [0,e], rj e = x/e — 1 on \e,2e\ and Xe = 1 on [2e. 1]. 

In a similar way to deriving (13.4.161) . we can get 


/ 


PeP 25 


2 V * + 2 


11 /l/\2 


dx + 4— j rje4>p 25 ^ l " > v 2 dx< 2 J p 2S Q 2 dx + 4 j P 2S dx, 


due to (I2.1.2p . Letting e -» 0 to give (|3.4.25p . Clearly, (I3.4.22p ~ A3.4.241) follow from similar argu¬ 
ments. 

□ 
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Lemma 3.15 Suppose that (|3.2.3D holds for a suitably small eo- Then for any a £ (0,1), 


Mli/x. < 2 \\v\\ ||u x ||, t£[0,T], 


\xv x \\ 2 L oo < c||xu x || (||x^ te || + ||u x || + \\v/x\\) , t £ [0,T], 


|||.°°([0,a]) - ( 1 / a ) II V x II L 2 ([0,a]) + 2 11 ^ II L 2 ([0,a]) \\ V xx\\L^{[Q,a\) > t € [°^]> 


(3.4.26) 

(3.4.27) 

(3.4.28) 


\ V / X \\l™([0,a]) < ( 1 / a )\\ V / X \\l^[0,a]) +2 ^ V / X ^L 2 i[0,a])^ V / X ^\\L^[0,a])^ ^ [°> T ]- (3.4.29) 


Proof. Clearly, (13.4.261) follows from the boundary condition v(0,t) = 0 and the Holder inequality. 
For xv x , notice that 


, \ 2 2 / V 'J 

[xv x ) = r x ( x— 


= 2r: 


fte) 


yj \ -y; y 


V — dy <c 


x- 


and 


= x(g t - 2-) <\\ x g tx \\ + 2 x(-) 

\r x J x V r/ x \rJ x 


Thus, 

H^llioc < c||am x || (||a;^te|| + ||v x || + \\v/x \\), 
which verifies (13.4.271) . ([3.4.281) and ([3.4.291) follow from simple calculations. 
Lemma 3.16 Let 5 be a fixed positive constant. Then for t £ [0, T], 

||r — x\\ 2 LOO < 2 II?’ — x|| ||r x — 1||, 

x 3/2 p 5 (r x - 1) <3 xp\r x - 1) +2 x 3 p 25 (r x - l)r. 


□ 




'2 

x 3 / 2 p 5 (r/x — 1) <3 xp s {r/x — 1) 


+ 2 


„3 - 25 , 



(3.4.30) 

^ L 1 

(3.4.31) 

- 1 ){r/x) x ^ . 

(3.4.32) 


Here (13.4.311) holds if the quantities appearing on its right-hand side are finite. 

Proof. Clearly, ([3.4.301) follows from ([2. 1.131) and the Holder inequality. For any x £ [0,1], 

x 3 p 25 (r x - l) 2 = 2 J ( y 3 p 25 (r y - 1 ?) y dy 

<3 I y 2 p 25 {r y - 1 ) 2 dy + 2 [ y 3 p 2S {r y - 1 )r yy dy, 

Jo Jo 

due to ([ 2.1. 2D . This gives (13.4.311) . Similarly, we can obtain ([3.4.321) . 


□ 


The following lemma is on the equivalence of the functionals £(t) and <£(t) and is the key to 
the verification of the a priori assumptions (13.2.31) and (I3.2.4D . 
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Lemma 3.17 Suppose that (|3.2.3D and (I3.2.4j) hold for suitably small numbers eo and e\. Then, 
\\(r x - l,r/x - l,v x ,v/x) (-,t)\\ 2 Lao < C£(t), t £ [0,T], (3.4.33) 

c£(t) < <S(i) < C£(t), t £ [0, T\. (3.4.34) 

Proof. The proof of ()3.4.33l) consists of two steps, in which the L°°-bounds on the intervals I\ = 
[0,1/2] and I 2 = [1/2,1] will be shown, respectively. Once ([3.4.331) is proved, (13.4.341) follows then 
from the definitions of £{t) and <£(t) by noticing that v(0,t) = r(0,t) = 0. 

Step 1 (away from the boundary). It follows from (13.4.31) . (I3.4.5P and (|3.2.1I) that 

\\Sxt\\ 2 < c (\\p^g x \\ 2 + \\pv t \\ 2 + || p{r - x,xr x - x)|| 2 ) < C£, (3.4.35) 

which, together with (I3.4.23P and (13.4.51) . implies 

\\p r ~ 1/2 (v xx ,(v/x) x )(-,t)\\ 2 < CS(t). 


Thus, 

\\( r xxi / X )xi v xxi Si C£(t). 

In view of (I3.2.9p . we then have 

||(r x - 1 ,r/x - l,v x ,v/x)(-,t)\\ 2 L 2 {h) < C\\(xr x -x,r- x,xv x ,v)(-,t )||| 2(Jl) 
+C\\(xr xx ,x(r/x) x ,xv xx ,x(v/x) x )(-,t)\\ 2 L3 ^ < C£(t). 

Hence, 

\\(r x - 1 ,r/x - l,v x ,v/x)(-,t)\\ 2 LOO ( h) 

<C\\{r x - 1 ,r/x - l,v x ,v/x)(-,t)\\ 2 H i( h) < C£(t). 

Step 2 (away from the origin). It follows from (13.2.ip and (I3.4.36p that 

II (r x - 1 ,r/x- l,v x ,v/x)(-,t)\\ 2 < C£(t), 

which implies, with the aid of (13.2.11) . that 

\\{r/x - l)(-,t)lli<*>(/ 2 ) < 411(7- ~ x)(-,t)\\ 2 Loo{h) < C\\(r - x) (-,t)\\ 2 H i (h) < C£(t ), 
||( 77 / 0 ;)(-, £)|| 2 oo(j 2 ) < 4||n(-, t)|| 2 oo(/ 2 ) < C \\v(-,t)\\ Hl (i 2 j <C£(t). 

Clearly, 

IK’"* - 1 )(-> i )lli~ (/2) < £(*)■ 

In view of (13.4.271) . (13.4.381) and ([3.4.351) . we see that 

IKM)ll£»(i a) < 4 II(^)(-,0IIl- ( 7 2) < 4\\( xv x)(-,t)\\ 2 L °° < C£{t). 

So, (13. 4. 331) is a consequence of (13. 4. 371) . (I3.4.39I) - ([3.4.42D . 


(3.4.36) 

(3.4.37) 

(3.4.38) 

(3.4.39) 

(3.4.40) 

(3.4.41) 

(3.4.42) 

□ 
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3.4.2 Part I: global existence and decay of strong solutions 

In this subsubsection, we prove the global existence and large time decay of the strong solution for 
suitably small <£(0). 


Lemma 3.18 Suppose that (13.2.31) and (13.2. 4p hold. Then there exist positive constants C, C(8 ) 
and C(8,a) independent oft such that for any 8 € (0, 2(y — l)/(37)) and a € (0,1), 


P 7 1/2 (r xx , (r/x) x ) (-,t) <C£( 0), 

(1 + t) (7 " 1)/7 " e || (r xx , (r/x) x ) (•,t)||i 2([M) < C(0,a ) (£(0) + || r 0x - l|||oc) , 

(1+ t) (7_1)/7_e p 1/2 v t (-,t) +[ (l + s) ( -' y ~ 1 ' >/ ' y - 0 \\(v x ,v/x,v sx ,v s /x)(-,s)\\ 2 ds 

Jo 

< C(8) (£(0) + ||roz — 1 ||l°°) > 

(1 + t)(7-i)/7-® ||(^, („/*).) (■,t)||i 2([M) < <7(0, a) [f ( 0 ) + ||r 0x - ini-] , 

/or all t £ [0, T]. 


(3.4.43) 

(3.4.44) 


(3.4.45) 

(3.4.46) 


Proof. The proof consists of four steps. 

Step 1. In this step, we prove (|3.4.43l) and ([3.4.441) . Multiplying equation ([3.4.3D by p 2 ^~ l Q x and 
integrating the product with respect to the spatial variable, one gets, using the Cauchy inequality, 
that 


dx. 


p2 ^ lgldx+ 1 / (?k) p3l ~ lglix 

<C J v 2 dx + C j x 2 p 7 — 1 j + (r x — l) 2 
It follows from (I3.4.47|) . ([3.3.121) and (|3.3.29l) that 

J p 27-1 f/ 2 (x, t)dx + J j p^^Gxdxds < C£( 0). 

This, together with (13.4.221) . gives ([3.4.431) . 

In a similar way to deriving ([3.4.471) . we can get 


(3.4.47) 


(3.4.48) 


lf t j^ eiix + 2j 

<C J v 2 t dx + C J x 2 p 1 


„2 \ 7 


r ^ r 


/ 7 2 g 2 dx 


- 1 ) + (fx — 1 ) 


dx. 


Multiply the inequality above by (1 + 1)^ 1 ^ 7 9 and integrate the product to give 

rt 


(i + t)7 9 J p 37 2 g 2 x (x,t)dx + 1 (1 + s)7 9 J j 


)^- 2 G 2 r dxds 


0 J (1 + s) t 9 J v 2 dxds + C J (1 + s) r 9 J x 2 p 1 -1 j +(r x — iy 

+ C J (1 + s)~J~ 9 Jp 3l ~ 2 Gldxds < C{8) (£(0) + ||r 0x — l|| 2 oo) , 


dxds 

(3.4.49) 
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due to (|3.3.41l) . (|3.4.48ll and the following estimate: 


J (1 + s) t 9 J p 37 2 Q^dxds 


= / (i+*r 


~^~ e 1 / a27-ir;2. 



p 21 1 Q%dx ) ( / p 37 1 Q 2 dx ) ds 


7~1 

7 


<C 


f (1 + s) 1 l9 ds sup f p 27 1 Q 2 (x,s)dx + f f p 3j 1 G 2 dxds. 

Jo sG[0,t] J Jo J 

It gives from ((3.4.2211 and (J3.4.49I) that 

(l + t)V " 0 f p 37 ” 2 (r 2 xx + \{r/x) x f) (x, t)dx < C(9) (S(0) + ||r 0x - l||£oo) 


(3.4.50) 


So, (13.4.441) follows directly from (13.4.501) . 
Step 2. In this step, we prove that 


J (l + s) 2 * 9 J (vl + \v/x\ 2 ^ dxds < C{0) (5(0) + \\r 0x - l|||oo) . 


(3.4.51) 


Multiplying the square of (13.4.31) by p 27 2 (1 + t)^ 7 1 ^ 7 6 and integrating the product with respect 
to the spatial and temporal variables, we have, using (|3.4.491) and (13.3.411) . that 

j (1 + s)~ 9 j p 2l ~ 2 Ql s dxds < C J (1 + s)~ 9 j (p 4l ~ 2 Gl + v 2 ) dxds 

+ C j (1 + s)~ 9 j x 2 p' i (^ - l) + (r x - l) 2 dxds < C(0) (5(0) + \\r 0x - l|||cc) . 

This, together with (I3.4.23P and (13.4.49p . implies 

J (1 + s) J S— 9 J(vl x + ^ dxds < C(0) (£(0) + ||ro* - l|||oc) . 


Clearly, 


^(l + s) 7 ^ 9 (vl x + ^ dxds <C(9) (£(0) + \\r 0x - l\\ 2 L oo) . (3.4.52) 


Then, it follows from (|3.2.9I) and (13.3.411) that 

rt .. , ri/2 

'0 Jo 


rt 7 -i r 1 / 2 

J o (1 + s)—~ e J [v 2 x + \v/x\ 2 j dxds < C{6) (£(0) + || r 0x - l||£oo) . 


We use (I3.3.4ip again to obtain (I3.4.51D . 

Step 3. In this step, we show that 

(1 + i)( 7_1 )/ 7_61 J pv 2 (x,t)dx + f(l + s)^- 9 f (y 2 s +{v s /x) 2 ) dxds 
<C{6) (£(0) + ||r 0a: -l||ioo). 


(3.4.53) 
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Differentiating (13.3.11) with respect to t yields 


_ ,'x\ 2 _ fx \ 3 v 

P - v tt -2p[-\ —v t — 7 


—C ( 2 — + — 

Vrr 




=P\— + 2 -) /' ( * 2 7o 


(3.4.54) 


Let ip be a non-increasing function defined on [0,1] satisfying 

ip = 1 on [0,1/4], ip = 0 on [1/2,1] and \ip'\ < 32. 

Multiplying equation (13.4.541) by ipvt and integrating the product with respect to the spatial vari¬ 
able, one has, using the integration by parts and the boundary condition w( 0 , t) = 0 (so vt{ 0 , t) = 0 ), 
that 


d r i 


dt / Y LI " 


X Y 2 


dx + p J ( ipv t ) x dx = J\ + J 2 4- J 3 , 


(3.4.55) 


where 


/ v { x \ 2 f /x \ 5 

-pip i—J v 2 dx + 4 / (—J vcppipvtdx < C / (n 2 + n/)dx, 

7 /(^) + 


J 2 :=- 

J 3 : =/ u 


u 2 n 2 
^2 ^71 

' X 1 


( ipvt) x dx. 


The second term on the left-hand side of (13.4.551) can be estimated as follows: 

{ipv t ) x dx = J ip^dx + 2 j ^v t v tx dx + J ip‘' 

V t dx ~ C J 0 ( x2 ^ + v t ) dx 


Vxt 2 Vt 
r x r 


Vxt | Yt 
r T r 


vtdx 


> [ i^dx - 
J r x 


> hP 


v v 2 


xt . T x Vj 
• o 
x V z 


r 1/2 

dx — C (x 2 n 2 t + n/) dx. 

Jo 


Then, 


d_ fl 

dt J 2 f 


(-) v^dx + p f ip 

' v xt , r x vl 

\rJ J 

r x r z 


r 1/2 

dx<C {x 2 v 2 xt + v 2 + v 2 ) dx + J 2 + J3. 
Jo 


It therefore follows from the Cauchy inequality that 


|/^(/) 2 »?*+f/v’ 


I o 


r 

This, together with (13.3.291) and ([3. 4. 511) . implies that 

rl/4 


/* 1/2 

dx<C { x 2 v xt + v 2 + {v/x ) 2 + v 2 ) dx. (3.4.56) 
Jo 


-7-1 P r 1 ^ rt 3-1 

(1 + i) t / pv 2 dx + (1 + s) 1 

Jo Jo 


-e 


V xs + “I ) (£(°) + ll r 0z - 1|||<X=) 
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Using (|3.3.29l) again, we obtain (13.4.531) . So, ()3.4.451) follows from (13.4.531) and (13.4.511) . 

Step 4 • I n this step, we prove (|3.4.46l) . It follows from (j3.4.3l) . (13.4.531) . (13. 3. 411) and (13.4.49p 
that 


p 7 2 Q 2 t (x, t)dx < C J pv 2 (x,t)dx + C j p 37 2 Q 2 {x,t)dx 


+ C / x z p 


2 tO 


+( r x-l ) 2 (x,t)dx < C(9)(l + t) 7 +e (5(0) + \\r 0x - l|||oo) 


This, together with (13.4.231) and (13.4.49^ . implies 

(1 + t)V" e y p 37 ~ 2 (vl x + \{v/x) x f} ( x,t)dx < C(9) (5(0) + ||ros - l||ioo) . (3.4.57) 

So, (I3.4.46P follows from (13.4.571) . 

□ 

Proof of Theorem 12.31 The proof of the first part of Theorem 12.31 follows from the local existence 
and uniqueness results of strong solutions, of which a sketched proof is given in Appendix, Part I. 
So, the global well-posedness of strong solutions with the estimate (|2.3.1I) can be shown by Lemma 
13.181 together with the equivalence of 5(f) and (£(t) shown in (|3.4.34D and the lower-order estimates 
obtained in Subsection 13.31 through the standard continuation argument . 

Moreover, (12.3.31) and (12.3.21) follow from (|3.4.100p and (13. 4. 1011) in Lemma 13.241 which will be 
proved later. 

□ 

To complete the proof of part i) of Theorem 13.11 it suffices to show the following Lemma. 


Lemma 3.19 For the global strong solution obtained in Theorem \2.3l there exist positive constants 
C(9) and C(9,a ) independent oft such that for any 9 £ (0, 2(7 — l)/(37)) and a € (0,1), 


(1 + t)(7-U/7-e || (ri _ i, r / x _ i, Vx ,v/x) (-,t)|| 2 < C(0)<£(O), (3.4.58) 

(1 + t ) (7_1)/7_0 11(77 - 1 ,r/x - l,v x ,v/x) (-, i )ll^i([o,a.]) ^ C(9,a)£( 0 ), (3.4.59) 

{1 + t) 2 h-i)h-0\\ {r -x)(.,t)\\U<C(9)m, ( 3 - 4 . 60 ) 

(l + t)^- 2 y^-^\(v,xv x )(-Mloo <C(9)m, (3.4.61) 

(1 + t) (7-i)/T-0 \\(v x , v /x) (.,f) llioo < C(0)<S(O), (3.4.62) 

(l + i)(7- 1 )/T-e ^(37-2)74^ _ l r / x _ i)( t \ 2 <C(6)£(0). (3.4.63) 

L°° 


Proof. It follows from (13.2.9p . (13.3.411) . (13.4.441) and (I3.4.46P that 

||(r7 - 1, r/x - l,T7j v/x) (-, t)||ia ([0jl/2]) < C\\{xr x - x,r - x, xv x , «)(■, t)||| 2([0)1/2]) 
+C\\(xr xx ,x(r/x) x ,xv xx ,x(v/x) x )(-,t)\\ 2 L 2 ^ 1/2 ^ < C(9)( 1 + f)” (7_1)/7+e (f(0), 

which, together with (13.3.4ip . gives (13.4.581) . Clearly, (|3.4.59p follows from (I3.4.58p . (|3.4.44D and 
([3.4.46p : (I3.4.60P from (I3.4.30p . (|3.3.4ip and (|3.4.58p : the estimate for v in (I3.4.6ip from (13.4.261) . 
(13.3.411) and (13.4.581) . Due to (|3.4.3D . (13. 3. 411) and (13.4.491) . we have 

\\xGxt(-,t)\\ 2 <C {\\xfg x {-,t)\\ 2 + \\xpv t (-,t)\\ 2 + \\p(r — x,xr x — x)(-,t)|| 2 ) 

<C(9)(1 + f) _ ^ 7_1 ^ 7+e C:(0), 
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which implies, using (13.4.271) . (13.3.411) and (13.4.581) . that 


(1 + t)(3?-2)/(2?)-* \\xv x (;t)\\l~ < C(d)£( 0). 

This verifies (13.4.611) . (13.4.621) follows from (13.4.611) . (I3.4.59D and the fact || • ||l°°([ 0 , 1 / 2 ]) < C || • 

Hff 1 ([o,i/2])- _ 

It follows from (13.3.41 ]) that 


xp h 1)/2 (r x - l)(-,t) < \\x(r x - 1 )(-, t )|| 1 7 

<C(0)(l + t)- (7 - 1) / 7+0/2 ^(O), 


xp l/2 (r x - 1 )(-,t) 


(7—i)/7 


which implies, using (13.4.311) . the Holder inequality and (13.4.431) . that 


x 3 / 2 ^( 3 7 2)/4(p x !)(•, t) 


< C 


xp { 7 1)/2 (r- a: - !)(-,/) 


P^ 1 1 )/ 2 r xx (-,t) 


+C || x(r x - 1)(-, t )|| 2 < C{9){ 1 + t)~^h +e m- 


This, together with (|3. 4.591) . gives 


-(37-21/4^ _ < C( 0)(1 +t)-('r- 1 )/ 7 +e (£ (o). 


L°° 


Similarly, we can use (13.4.321) to get 


rf?m-W(r/ x _ i)(.,t) < C(0)(1 +i)-(7- 1 )/7+e (£ (o). 


L°° 


This finishes the proof of (13.4.631) . 


□ 


3.4.3 Part II: faster decay 

In this subsection, we prove part ii ) of Theorem 13.11 under the assumption 

&«(0) < 00 , a € ( 0 , 7 ). (3.4.64) 

The estimates in this subsection are for the global strong solution of (12.1. Ill) as stated in Theorem 

12731 

To obtain the faster time decay estimates of the higher-order norms, we rewrite equation (13. 4. 3D 
in the form of 


x 


fP(x, t) + nQxt = —pv t , where <P(x, t) : = 7 


2 - \ 7 

x z p X 


r 2 r x 


Gx + 


„2 \ 7 


r 2 r x 


- 1 -) 4 

. r 




It should be noted that 

x‘ 


% = IP 1 


r 2 r 1 


Gxt+ty 1 , where := 7 p 7 


x 


2 \ 7 






X 


2 \ 7 


r 2 r 


- 1 7 ) 


'*/ J t 

This equation is convenient for us to derive the time decay estimates for r xx and v : 


x(f>p. (3.4.65) 

x<j>p. (3.4.66) 
with weights. 
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Lemma 3.20 Let a £ (0,7) and 3a(0) < 00. For the global strong solution obtained in Theorem 
2.3 1 there exist positive constants C(a ) and C(9,a) independent of t such that for any 0 < 0 < 
min{2(7 - l)/(37), 2(7 - oO/y}, 


^(27 1 a )/2 ( r ^ (r/x) x ) (-,f) < C'(a)3 r Q (0), 

(1 + t )(7-i)/7-« p(^-)/\ rx _ 1, r /x - 1)(-, t) 2 < C(0, a)ff Q (0). 

L°° 

Proof. In a similar way to deriving (I3.4.48p . we have 

J p 2l ~ 1 ~ a C} 2 (x 1 t)dx + J j p 3l ~ l ~ a Q 2 dxds < C j J v 2 dxds 


(3.4.67) 

(3.4.68) 


+ C 


II 


(3.4.69) 


x 2 p l+l - a 


( r/x - l) 2 + (r x - 1)“ dxds < C'(a)3' a (0), 


which, together with (13.4.221) . gives (|3.4.67|) . The proof for (13.4.681) is the same as that for (13. 4. 631) . 
We omit the detail here. 

□ 

Lemma 3.21 Let a £ [7 — 1,7) and 3»(0) < 00. For the global strong solution obtained in 
Theorem 12.31 there exist positive constants C(a,9) and C(a, 8 ,a) independent oft such that for 
any 0 < 8 < min{2(7 — 1)/(37), 2(7 — a)/ 7} and a £ (0,1), 


K I 47 — 3 3 A 

o o-,, o u 


(1 + f) 2+ 2 


p 1/2 v t (-,t) +[ (1 + s ) 2+2 7 2 ° \\(v x ,v/x,v sx ,v s /x) (-,s)\\ 2 ds 

Jo 

<C(a,8)(£( 0 ) + l)&,( 0 ), (3.4.70) 

min ( K I 4 ^~ 3 3 P k. j IO-y- 9 9nl „ 

(! + *) 27 2 ’ 4 47 4 ^ \\{rxx,('r/x) x ,v xx ,(v/x) x )(-,t)\\ 2 L 2 {[0a]) 

<C(a,9,a)(£(0) + l)$ a (0). (3.4.71) 

Here k = 0 when a = 7 — 1, and k = ( 1 / 7 ) min{a — (7 — 1), 7 — 1} — 8 when a £ (7 — 1 , 7 ). 

Proof. The proof consists of four steps. 

Step 1. In this step, we prove 


(1 + t) K J p y Ql(x,t)dx + J (1 + s) K J p 2l Gldxds < C( 8 , a)3h(0), 
J (1 + s) K j (Q 2 x + v 2 + \v/x\ 2 ) dxds < 0(8, a)3a(0). 


(3.4.72) 

(3.4.73) 


When a = 7 — 1, ([3.4.691) implies (13.4.721) . When a £ (7 — 1,7), we use the same way as that 
for the derivation of ([3.4.491) to obtain 

(1 + *)-/ P 1 Qx( x -,t)dx + J (! + ,)«/ p 2l Qldxds 

(1 + s) K J v 2 dxds + cj (1 + s) K J x 2 p 1 - l) +(r x -iy 
+ C J\ 1 + s)*" 1 J fPGldxds < C(8, a)&,(0), 


<C 


dxds 
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due to 


and 


\ (27-l-a)/7 / r \ (a+l-7)/7 

d 7 ^x < ( / p 27_1_Q ^dxJ (J p^^ l ~ a Qldx J 


J (l + s) K 1 J p^G 2 x dxds<C J I p 37 1 a Q 2 dxds 

+ C [\l + s) -^-^ds sup / p 27_a £/ 2 (x, s)dx. 
J 0 sG[ 0 ,t] J 

This verifies ()3.4.721) . It follows from (13.4.31) . (13.4.721) and (J3.3.41D that 

J (1 + s) K J Q 2 x dxds < C(0, a)5a(0), 


which, together with (13.4.51 p . gives (I3.4.73p . 

Step 2. In this step, we prove that 

(l + t)2 + 2 t ~2 e J x 2 p^- 2 g^{x,t)dx + J (l + s)2 +_ 2r z e J x 2 p 3l ~ 2 Ql s dxds 

< C(a,Qy$ a (0), where &*(0) := 3a(0) + £(0)£ Q (0). 

Multiplying (|3.4.65D by x 2 p 27 " 2 < }3f and integrating the resulting equation give that 

J x 2 p 27 ~ 2 fp 2 dx + p J x 2 fj 2l ~ 2 %Gxtdx = J x 2 p 2 'i- 1 y t ^v t dx, 

which implies, using (I3.4.66p . that 

1 ^ ^ c 2 p 27 _ 2 < p 2 dx + p 7 J x 2 p 3l ~ 2 


1 d_ 

2 dt 


2 dt 


„2 \ 7 


r 2 r I Gxtdx 


=7 f x 2 p 37 1 Gxt%v t dx+ f x 2 p 21 1( Pi ^v t dx-p (x 2 p 21 2 ^iG xt dx. 

J \r z r x J r z J r z J 

Thus, one has 

\-JiJ x2 p 2l ~ 2 ^ 2dx +I ~y)' x 2 P 3i ~ 2 (^t) <C j (x 2 p 7 ~ 2 <Pf + u t 2 ) dx 


< C J ( x 2 v 2 + v 2 ) p 37 2 G 2 dx + C J (x 2 u 2 + v 2 + u 2 ) dx. 
Combining this with (13.4.491) shows that 


i 4 / xV -V^ + 

<C fHxiUlr 


?/ 


*V 7 - 2 , „ 2 


„2 \ 7 


r 7 


QUx 


[\\xv x \\ 2 Laa + C \\v\\ 2 La< ^j J p 37 2 Q 2 dx + C J (x 2 vl + v 2 + v 2 ) dx 
<C(0)<£(0)(1 + i)~(7-!)/7+f ^Hxu^ll^oc + ||uH^oo) + Cy (x 2 u 2 + u 2 + u t 2 ) dx. 
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(3.4.74) 


(3.4.75) 
























It then follows from (13.4.271) and (13.4.261) that 


i d r 

2dt J 


x 2 p 27 2 ip 2 dx + 


?/ 


x 2 p^~ 2 


'2 \ 7 


r ^ r 


<C(0)g(O)(l+t)^-^ + ^ / (C + ^ + W*) 2 )dx 


+ C(6 , )G:(0)(1 + t) 2 ^ 2 + 2 J {x 2 v 2 + v 2 )dx + C J (x 2 v 2 + v 2 + v 2 ) dx. 

K I 47 — 3 3 g 

Multiplying the inequality above by (1 + i) 2 2 r 2 and integrating the product give that 


K, 1 47 — 3 3 /) 

(1 +f) 2+ 2 

<c 


' j x 2 p 2,1 2 ?($ 2 (x,t)dx + J (l + s) 2+ 2 t .■/ x 2 /? 37 2 Q 2 s dxds 


J x 2 p 27 2 fp 2 (x,0)dx + C(d)<£(0) J (1 + s) K J (Q 2 x + x 2 + (x/x) 2 ) dxds 


K 1 47 — 3 3/T 

2 ' 27 2 17 


+ C 

Note that 


+ J [c(0)<£(O)(l + s)V- fl +(7(1 + 8) 

j\l + S )d + ,)- +2 ^- e J x ip>i-^ dxds . 


x 2 v 2 + v 2 + v 2 ) dxds 


(3.4.76) 


J x 2 p 2 ' 7 2 ip 2 dx < C J x 2 p Al 2 G 2 dx + C J x 2 (P — 1^ + (r x — l)"' 

* + flzi _ < 23^1 _ e 

2 27 2 — 7 ’2 27 2-7 


dx, 


(3.4.77) 


Then, it follows from (13.4.761) . (|3.4.73p . (I3.3.41P and (|3.4.49l) that 

(l + t) 2 ^ 27 20 J x 2 p 2l ~ 2 ^ 2 dx + J (l + s) 2_l 27 2 61 J x 2 p 3l ~ 2 Gl s dxds < C(a,9)$ a ( 0). 


This, together with (13.3.411) . implies (13.4.741) . 
Step 3. In this step, we prove that 


J (l + s) 2 + 2 ¥ 20 J (v 2 + (v/x) 2 ) dxds < C(a,0)$ a (0 ), (3.4.78) 

w, 1 47 — 3 3/j f n ft w, 1 47 — 3 3/1 /" , ,, n 

(1 + t) 2 2 t 2 / pv 2 (x,t)dx+ / (1 + s) 2 27 2 / (+ 2 S + (x,;/x) 2 ) dxds 

<C(a,(9)5a(0). (3.4.79) 


As a consequence of (13.4.781) and (13.4.791) . we get (13.4.701) . 
It follows from (13.4.lip . (13.2.91) and (13.4.61) that 


f (v 2 + (v/x) 2 ) dx < 4 [ g 2 dx + 4 [ G 2 dx<c f x 2 (G 2 + Q 2 X ) dx + 16 f x 2 G 2 dx 

J ‘ JO dl/2 do d 1/2 

< C J x 2 G 2 dx + C J x 2 p 3 ' y ~ 2 Gt x dx < C J (x 2 u 2 + x 2 ) + C J x 2 p 3l ~ 2 G 2 x dx- (3.4.80) 
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This, together with (13.4.741) . (|3.4.77|) and (13. 3. 411) . gives (13. 4. 781) . With the aid of (J3.4.78D . we can 
use the same way as to the derivation of (13.4.531) to obtain (j3.4.79l) . 

Step 4■ In this step, we prove 


• / k i 47 — 3 3 n k, 1 IO 7 —9 9 n\ /* ~ 

(l + t) mm l 2 + 2 7 2 e - 4+ 4 7 4 '/ p4 7 -2^2 fat)dx < C{a,6)3 a (0), 

min f K I 4t ~ 3 3 P K I 10 7~9 9q1 I" ~ 

(l + t ) mm i 2+ 2 7 2 e ’ 4+ 4 7 < C(a,0)$ a (0). 


With (13. 4. 811) and (|3.4.82D . we can obtain (13.4.711) by use of (13.4.221) and (13.4.231) . 
Let be a non-increasing function defined on [0,1] satisfying 

= 1 on [0,1/8], -0 = 0 on [1/4,1] and |0| < 32. 

Following the derivation of (13.4.751) . one can obtain 

# 37 " 2 






<C , (0)C(O)(1 + t) (7 1)/7+0 ||(^,'y/^)|lioc ([0il/4 ] ) +cj (v 2 x + {v/xf + V?) dx. 


(3.4.81) 

(3.4.82) 


In view of (13.4.281) and (13.4.291) . we see 

II II 2 L II / II 2 

IIW|Il°°([0,1/4]) + \\ v l x \\L°°{[0,\/i]) 

<C 


r , . r r / N 1 V 2 ri / 2 / \ 

J (v 2 x + \v/x\ 2 )dx + C j (v 2 x + \v/x\ 2 ^dx JI 


which implies, 

5 s/ + f / ^- 2 (^r) 

5-6 k 5 7 1 / 2 

<C(0)<£(O)(1 + t)^7 I_4+20 / (uL + Kv/®)*! 2 ) dx 

■Jo 

+ C'( 0 )(£(O )(1 + i) - ^ J ( v x + ( v / x ) 2 ) dx + C J (v 2 + (v/x) 2 + v 2 ) dx. 

Similar to (13.4.741) . one can use (|3.4.78D and (13. 4. 521) to obtain 


k 1 10 7 —9 9/1 f — . 00 0 re 

(1 + t) 4 4 t 4 / 'ipp 1 ~ 2 Q 2 (x,t)dx + / (1 + s) 4 

<C(a, 0 )(L(O), 


10 7 —9 9p 
4 7 4 


-37-2^2 


Q xs dxds 


due to 


k IO 7 — 9 

T ^- a - 

4 4y 



27-I 

7 


0 . 


This, together with (I3.4.74p . implies (I3.4.8ip . Finally, we can use (13.4.31) . (I3.4.79p . (I3.3.41P and 
(13.4.8ip to show (|3.4.82l) . 

□ 
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Lemma 3.22 Let a £ [7 — 1 , 7 ) and 3a(0) < Tor the global strong solution obtained in 
Theorem, 1 2. ‘A there exist positive constants C(a,9) and C(a,9,a) independent oft such that for 
any 0 < 0 < min{ 2(7 — 1 )/( 37 ), 2(7 — a)/^} and a £ ( 0 , 1 ), 


a + ty- 


K 1 47-3 3 


27 


-%e 


(r x ~ l,r/x - 1) (■,i)||i 2 ([0iO]) < C{a,9,a)$ a (0), 


(1 + t) 2+ 27 2 e \\{v x ,v/x) (',t )|| 2 < C(a,9)3 a (0), 

min f K I 4 T'~ 3 3 p k | 107-9 9/0 „ 

(! + *) 27 2 ’ 4 47 4 J \\(r x -l,r/x-l,v x ,v/x)(-,t)\\^ M) 

< C(a,d,a,y$ a (0), 

(1 + t)^+f-! 0 ||w(.,i)||| ao < C(a,0)^(0), 

(1 + t)^ bl \\xv x (-,t)\\ 2 Lao + (1 + t)? mm{bl ’ b2 -' \\{y x ,v/x) (•,i)|| 2 oo < C{a,0)$ a {G), 


H. _i_ 27—1 _ §_ 

(1 +t) 2+ 2 7 2 


f /2 {r x - 1, r/x - 1)(-, t) < C(a, 0)$ a (O). 

L°° 


(3.4.83) 

(3.4.84) 

(3.4.85) 

(3.4.86) 

(3.4.87) 

(3.4.88) 


Here 3^(0) = 3a (0) + £(0)3a(0), k = 0 when a = 7 — n = ( 1 / 7 ) min{a — (7 — 1), 7 — 1} — 9 
when a € (7 — 1 , 7 ), 


l K 47 — 3 3 

j\ = mm < max <1 —|---— 9 


2 7 

1 _i_ 4 7 ~ 3 _ 3 0 

2 + 27 2 


+ 


2 

27-I 

7 


a + 1 3 , 2 7 -l 9 

27 - 1 + a’ 2 K + ~ 2 ^ 2 




i.2= m in{-+ ^ 


k 47 — 3 3 k IO 7 — 9 9 


2 7 4 47 

Proof. It follows from (13.2.9P and (|3.4.4|) that 

-1/2 

+ 4 

/o 4 1/2 

<c r " 


-7« f-f 


K 47 — 3 3 


27 


[ p y G 2 dx <C f Q 2 dx + 4 / x 2 p" l G 2 dx <cf x 2 [C, 2 + £ 2 ) dx + 4 f x 2 p y G x 
J Jo J 1/2 io 41/2 

/ x 2 (PQ 2 dx + C J x 2 p Al ^ 2 G 2 dx 
<C j ff 1 {\xr x — x\ 2 + |r — x| 2 ) dx + C J x 2 p Al ~ 2 G 2 dx, 


(3.4.89) 

(3.4.90) 

dx 


which implies, using ([3.4.241) . (13. 3. 411) and (13. 4. 741) . that 

(1 + t) 2 J p 7 (|r x - 1| 2 + | r/x - 1| 2 ) (x, t)dx < C(9, a)3 a (0). 
This gives (13.4.831) . It follows from (13.4.31) . (13.3.411) and (13.4.741) that 


(l + t f 2 +^-¥ j X 2 p 2ry - 2 g 2 x (x,t)dx < C(a,9)da( 0 ), 

which, together with (13. 4. 801) and (|3.3.41D . gives (I3.4.84|) . ([3.4.85P is a conclusion of (13.4.831) . 

(13.4.841) and (13.4.711) . (13.4.861) follows from (13.4.261) . (13.4.841) and (13.3.411) . 

Next, we prove (13.4.871) . It follows from (13.4.741) and (13.4.691) that 


(1 + 1 ) 


47-3 3 n\ Q + l 

27 2 U ) 27-I +ot 


J x 2 p 2l Ql{x 1 t)dx < C(a, 9)$ a (0), 


(3.4.91) 
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due to 


/■ 


Q + l 

27—1 + CK 


x 2-2y-\-ag2dx 


27-2 

27 —l+o: 


x 2 p 2 ' y Gldx — yj ^ 2 p 47 2 <5 2 d£ 

In a similar way to deriving (13.4.741) . we have 

(1 + t) 5K+ 2 t ~ 2 J x 2 p 2l G 2 (x,t)dx + J (1 + s) 5KH *7 2 J x 2 p^G 2 s dxds < (7(a!, #)5a(0), 
since (13.4.751) can be replaced by 


1 d [ 2 m 2 , , M7 /' 

--j xVdx+T j 


M7 / 2-7 
a; p' 


t 2 \7 

^ 2 t dx 


<<7 llsn 


■'aMlL 0 


+ < 


'C|Mli°°) j fG 2 dx + C J (x 2 v 2 + v 2 + v 2 )dx 

<C'(0)2(O)(1 + f) _K ^lIxu^H^oo + Hull^oo^ + C J {x 2 v 2 + v 2 + v 2 ) dx. 

This, together with (13.4.911) . gives 

(1 + t) max {(t + ^r-f 0 )2y^- J x 2p^g*( Xj t)dx <C(a,0y5 a (0). 

So, it yields from (|3.4.3p . (13.4.701) . (|3.4.92l) and (13.3.411) that 

(1 + t) bl — +9 Jx 2 Gtx( x ,t)dx < C(a,9)$ a (0), 

where is defined by (|3.4.89l) . This, together with (13.4.271) . (|3.3.4ip and (13. 4. 841) . gives 

(1 + t) bl/2 \\xv x (-,t)\\ 2 Lao < C(a,9yg a (0). 

Moreover, it follows from (13.4.281) . (13. 4. 291) . (13.4.841) and (I3.4.85P that 

{1 + t) b2/2 \\(v x ,v/x)(-,t)\\ 2 Lxm/4]) < C(a,9)$ a (0), 

where 62 is defined by (|3.4.90l) . Therefore, we have (I3.4.87p . 

In a similar way to deriving (I3.4.94p . one can get 

(1 + t) b2/2 \\(r x - 1, r/x - 1) (■) t) 11 Z/Oo ([0 ,i/ 4 ]) < C(a, ^)5a(0), 

where 62 is defined in (13.4.901) . It follows from (13.4.311) . (13.3.411) . (13.4.221) and (I3.4.72P that 


(3.4.92) 


(3.4.93) 


(3.4.94) 


(3.4.95) 


s 3 / V / 2 (r x -l)(-,t) 


< C 


xp l/2 (r x -l)(-,t) p 7/2 r ra (-,t) 


+C 


xp l/2 (r x - < C(a,9)$ a (0)(l +t) 2 27 + 2 . 


Clearly, we can get the same estimate for r/x — 1. Due to k + (2y — 1 )/ r ) — 9 < 62 , we then obtain 
(13.4.881) . 

□ 


Remark 3.23 As a consequence of Lemma 17. 22[ we have for any 0 <9 < 2(7 — l)/( 37 ), 

117-10 5 n O ~ 

(1 + f) 4 7 2 \\{v xx ,{v/x) x ) (•,t)lli 2 ([0 ,i/ 2 ]) < C(9)^- 1 e(0), 

(1 +t)^ 20 \\(XV X ,V) (-,t)|| 2 oc < <7(19)57-70(0), 


57 — 4 


(1 +t)-^~ 26 II xGtxMf < C(0)5 7 - 7 e( 0). 


(3.4.96) 

(3.4.97) 

(3.4.98) 


Indeed, (13. 4. 961) follows from (13.4.851) . (|3.4.97f) from (|3.4.86l) and (13.4.871) . and (13.4.981) from 
(13X931) . 
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3.4.4 Part III: further regularity 


In this subsection, we further study the higher regularity of the strong solution obtained in Theorem 
12.31 and prove part in) and part iv) of Theorem 13.11 

Lemma 3.24 Let a £ [ 7 , 27 — 1] and 5a(0) < 00 . For the global strong solution obtained in 
Theorem \2.3l there exist positive constants C and C{6) independent oft such that for any 0 <9 < 
2(7 - l)/( 37 ), 


-(2 7 -l-a)/2 (r ^ {r/x)x) ( . >t) < C £ q(q) + C{0){1 + i) (a- 7 +fl7)/(a-l) (£(0)j 


Ij_ 


{fxxi (r/x) x ) < ^ 2 7 _i( 0 ) + <7(0X1 + t)-^—*°£( 0 ), 


(3.4.99) 

(3.4.100) 


_7t-6 


IIK*, (v/x) x ) (-,t)f < C(9)( 1 + t)-T^+ 4 % 7 _ 1 ( 0 ) (1 + J 27 _i(0)) (1 + <£(0)), (3.4.101) 

provided that §' 27-1 (0) <00 in (13.4.1001) and ()3.4.101l) . 

Proof. In a similar way to the derivation of (|3.4.48D . we have 


f^- l - a gl( x ,t) dx + J* J f'<- x - a Qld.,xds 
* a (0) + C'jT J 


xV +1 ~“ 


(r/x - if + (r x - if 


dxds. 


It follows from the Holder inequality that 


/ 


< 


-7+1-a 


(r — xf + ( xr x — xf 


(r — xf + ( xr x — xf 


dx 

(a- 1)/ 7 


dx 


p 1 


(r — xf + (xr x — xf 


< (1 + t)-^~ l)h+e ^(1 + t)^- l)h ~ e J [(r - xf + ( xr x - xf 

/ r r -1 \ (7+1 -a)h 

x ( (1 + t)h>' -1 )/T' -£l / fp [r — xf + ( xr x — xf dx j 

which, together with (13.3.411) and the Young inequality, implies that 


dx 


dx ) 
(a-l)/ 7 


\ (7+l-«)/7 




a 


X 2 p< +1 - a 


(r/x - 1)" + (r x - 1)" 


dxds 


< C 


/w 

Jo 


7 -l)/ 7 - 6 » / - 7 


(r — x) 2 + ( xr x — x) 2 


dxds 


+ C(1 + t)i-(7-i-7«)/(«-i) sup j(i + a )(7-i)/7-e f [(r - xf + ( 

se[o,t] l J 

< C , (0)g(O)(l + t) (a_7+0 ' Y)/(Q “ 1) . 


X7V — X 


(x,s)dx 


Thus, 


J p 2 ^- 1 - a gl(x,t)dx + J* J p^- 1 - a gldxds<C 3 a {o) + c{d)£{o)(i + t)' L ^ x . ( 3 . 4 . 102 ) 
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This, together with (13.4.221) and (|3.4.13l) . gives (13.4.991) . 

Choose a = 27 — 1 in (13.4.1021) to give 

J g 2 x {x,t)dx < Cfo 7 -i(0) +C(0)(l + i)^V(O). (3.4.103) 

This, together with (13. 4. 131) . gives (13. 4. 1001) . It follows from (13.4.211) . (13.4.1031) . (J3.4.97I) and (13.4.981) 
that 

||( u xxj ( V / X )x ) (') t) Hl 2 ([i/2,1]) 

< c { 6 %_ ie { 0) {(1 + ty^ +2B + (1 + (fo 7 _1(0) + (1 + ty + ^ e(o))} 

< C ( 0 %_ ie ( 0) (^27-1(0) +1) (1 + tr^ +4e . 

This, together with (13.4.961) . implies (13.4.1011) . 

□ 

Lemma 3.25 Suppose that \\xp l ^ 2 vtt{-, 0) || 2 < 00. For the global strong solution obtained in The¬ 
orem 1 2. ,'A there exists a positive constant C independent of t such that 

poo 

\\xp 1/2 v u (-,t)W 2 + / \\(v ss ,xv ss x)(-,s)\\ 2 ds <C£(0)+ C\\xp 1/2 v t t(-,0)\\ 2 , t>0. (3.4.104) 

Jo 

Proof. Multiplying <9? (13.3.11) by r 2 vtt and integrating the resulting equation both in x and t. one 
can show (13.4.1041) . Indeed, the derivation of (|3.4.104l) is similar to that of (13.3.381) . so we omit the 
details here. 

□ 


4 Proof of Theorem 12.4 


Due to Theorem 13.II and Theorem 12.31 the triple ( p,u,R(t )) (t > 0) defined by (|2.3.41) and (|2.3.5I) 
gives the unique global strong solution to the free boundary problem (|1.1.2D . The decay estimates 
i) and ii ) in Theorem [2~4l follow from the corresponding ones in Theorem 13.11 by noting that 

| p(r(x,t),t) - p(x )| < Cp(x) (| r x (x,t) - 1| + \x~ 1 r(x,t ) - 1|) 


and 

v x (x,t) u(r,t ) x v(x,t) 

Ur{r,t) = and -= —tt-• 

r x (x,t) r r[x, t) x 

The Ty 2,00 -estimate of Rft) can be proved as follows. First, it follows from (13.3.291) that 

poo poo pi 

J v 2 (l,t)dt<C J J^ (v 2 + v 2 t )dxdt < C<£(0). 

One the other hand, (13.1.101) implies that 

poo poo pi 

/ «tt(M)di / (4 + vltt)dxdt < C£{0) + C\\xp 1/2 v tt (-,0)\\ 2 . 

Jo Jo Vi 
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Combining these two estimates with the fact that 


R\t) 


VfiX, t) < 1^(1,0) + 2 



/ roo \ l/ 2 

U v tt( l -.t)dtj 


gives (12.3.161) immediately. This finishes the proof of Theorem 12.41 


□ 


Remark 4.1 One may prove the boundedness of ru{x,t) for any fixed x G [0,1] and t > 0 if 
|vt(x,0)| is finite by an argument similar to the above. This implies that every particle moving with 
the fluid has the bounded acceleration for t G (0, oo) if it does so initially. 
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Appendix 


Part I. Local existence of strong solutions in the functional space (£(£) 

In this part, we prove the local existence of strong solutions to problem (12.1.11ft on a time interval 
[0,T*] for some T* > 0 in the function space {(r,v) : (£ £ C*([0, T*])} by using a finite difference 
method, as used in [HU [251 ESI E|, where either a one-dimensional model or a three-dimensional 
model with spherical symmetry in a cut-off domain excluding a neighborhood of the origin is 
considered. Ideas used to derive the estimates in Section [3] will be employed here to deal with the 
differences between the problem considered in this paper and those considered in [321 [23 [261 E]- 
The proof works for the case when ||r' a ,(ic, 0) — 1 ||l°° is small. It may be possible to obtain a 
local existence theory in our function space by only assuming that r x (x, 0) has both lower and 
upper positive bounds. However, we do not pursue this generality here which may need extra work 
because the main purpose of this paper is to establish the global existence of strong solutions for 
small data. 

Recall that (12.1.llaft reads 

K~r )’”‘-(® + 4A u) I = {^ 

where 53 is given by (12.1.121) . that is, 



We use (r°,v°)(x) to denote the initial data (r, v)(x, 0)(This notation will avoid possible confu¬ 
sions when we define the finite difference scheme.) The finite difference scheme is defined as follows. 
Let N be a positive integer, and h = 1/N. For n = 0,1, ■ • • N, set 


i -1 ^4 

tv 


+ ~ 1 


x 


(A-l) 


x n = nh, p n = p(x n ), q n = ( p 1 ) x {x n ). 

We can approximate (IA-11) by the following initial value problem for the system for v \, • • • 
of ordinary differential equations: 


•En 
Pn — 


dv n 1 
dt h 


53 n + 4A 


Di—1 

1- 

r n -1 


- ( 53 n +i + 4Ai — 




+ 


1 


Pn+l 


1 - 


k K 

Dx(0) = v°(x n ), 


1~n+l 1"n 


7 /„. \ 2 7 
n 


pi 


1 - 


r n ~ r n -1 


n , 

7 / T , \ 27' 

^n —1 

r n -1 


(A-2) 
, VN-l 


(A-3a) 

(A-3b) 
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where 


r n {t) = r°(x n ) + / v n (s)ds , 
Jo 


®n = A* 


^n —1 

r n - r n -1 


+ (2A2 — 77A] 


Dt-l 

1 


(A-4) 

(A-5) 


This system is supplemented by the following conditions to match the boundary conditions v(0,t) = 
0 and 35(1, t) = 0: 


vo{t) = 0, 35 N (t) = // 


VN ~ VN-1 
r N - r N -1 


+ ( 2 A 2 — —A 


VN-l 

r N—l 


= 0 . 


(A-6) 


Clearly, if r(0, 0) = 0, it follows from vq (t) = 0 that 

r 0 (t) = 0. (A-7) 

Next, we use the condition 35 at(A) = 0 to determine vn and r jv in terms of ujv_i and rjv- 1 . It 
follows from ( d/dt)r n = v n that 


0 = 35jv(t) = ^ { At In 


r N ~r N - 1 
h 


+ ( 2A2 — -Ai ) In 


Ov-i 

X N -1 


which implies 

rN{t) = r N -\(t) + [r°(x N ) - r°(x jv— 1 )] 
This, together with (1A-6I) . gives 


2A 2 

r°(x]y-i)\ M 3 M 

J'jv-i(i) 


uat(A) = v N -i{t) - f— - [r°(x N ) - r°(a:Ar_i)] . 

V At 3 p J 1 V Dv-i(o 

With (|A-6D - (IA-9|) . ODE system (IA-3D is closed. 

The approximation of the functional (£ dehned in (|2.2.6I) is given by 


2A 9 4Ai 


r°(z7V-l)^ M ^jV-l(t) 

DV-l(t) ' 


(A- 8 ) 


(A-9) 


&N(t) = max 
l<n<N 


r n {t ) - r n _i(t) 2 v n (t) - v n —\(t) 
h + h 


N—l 


+ h ^ Pn\Vnit)^ 


n=l 


N—l 


+ h^2 


rn 


n= 1 


r n+ i(t) - 2 r n (t) + r n -i(t) 


K 2 


+ 


1 / r n (t) r n _i (A) 


h \ x r 


%n—1 


(A-10) 


where and thereafter, <? = ( d/dt)g for a function 5 = g(t). It should be noted that h Yln=i Pn\v n (0)\ 2 
can be given by the initial data via equation (IA-3al) . From now on, we choose a positive integer Nq 
so large that 


(£jy(0) < 2<£(0), N > Nq. (A-ll) 

Lemma A-2 Let 7 £ (4/3,2). Suppose that <£(0) < 00 and r°(0) = 0. Then there exist positive 
constants N\ > Nq and e > 0 independent of N such that problem (1A-3D admits a unique solution 
(fn,Vn)(t) on [0, T*] for some positive constant T* independent of N satisfying 

<Sjv(t) < K£ n ( 0) < 2K<E(0), t £ [0, T*], N > N\ (A-12) 
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for some positive constant K independent of N, provided that 


k2-i| 




Moreover, T* satisfies 

rr* ^ ■ S c 1 

\ y/2K£(0) ’ K(1 + 2K£(0)) 

for some positive constant c independent of N. 


(A-13) 


Remark A-3 7 G (4/3,2) is in general not necessary for the local existence, 7 > 1 should be 
sufficient. The reason we put this condition in the lemma is to ensure the existence and uniqueness 
of the stationary solution, the Lane-Emden solution, to keep the consistence with the global existence 
theory. 


Proof of Lemma IA-SI It follows from the ODE theory that problem (IA-31) has a solution on a time 
interval. Let T/v > 0 be the maximum existence time. It follows from vq (t) = 0 that 


max \v n (t)/x n \ < £ N (t), t € [0,T N )\ 

l<n<N 


(A-14) 


which, together with r n = v n , implies that for n = 1, • • ■ , N, and t G [0, T/v), 


r n {t) - r n -i(t) 


h 


- 1 


< 


rp(x n ) - r 0 (a; n -i) _ ^ F* 

h Jo 


V n (s) - fn-l(s) 


h 


ds < ||r° - l|| Loo(J) +t sup a/^a r(s). 


se[o,t] 


This implies, due to ro(t) = 0, that for n = 1, • • ■ , N, and t G [0, T/v), 


\r n (t)/x n - 1| < ||r° - l\\ Loo(I) < || r° x - l\\ x(I) +t sup y/<£ N (s). 

w w sG[0,t] 


Therefore, one can check that, for n = 1, • • • , N, and t G [0, T], 
\(r n (t) - r n _i(i)) /h - 1| < 2e and \r n (t)/x n - 1| < 2e. 
provided 


T sup y/ <£n(s) < £, 
s€[0,T] 


and 


I r ^ — ill ^ f 

V x -*-11 Z,oo(/) ^ fc 


for some constant e > 0. In particular, it holds that for n = 1, • • • , N, and t G [0, T], 


1 r n (t) 3 1 r n (t) — r n -i(t) 3 

_ < _n W < _ and _ < -n_LW < _ 

2 ~ Xn ~ 2 2 - h “2’ 


if 


1 

6 ~ 4' 


(A-15) 

(A-16) 

(A-17) 

(A-18) 

(A-19) 
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With (1A-141) . (IA-151) and (IA-181) . we will prove in Steps 1 -3 that for sufficiently large N, there exists 
constant C > 1 independent of N such that if T satisfies (IA-16D . then 


01a+) 8 (70:+0) + C f [ 0 :tv(s) + *^Af( s )] ds, t £ [0, T]. 
Jo 


(A-20) 


Once this statement is proved, the lemma will follow from an argument which we give in Step 4- 
Step 1. In this step, we prove that 


TV—1 


r t N-l 


n= 1 


n =1 


( h n+ i(s) - v n (s) 


h Y + f hY^ X 2 n ( ' " + U °^ U " A *' ) + Vn(s) 

Jo 




h 


ds 


<C& N (0) + C [ [ 0 ;+s) + ^at( s )] ds. 
Jo 


(A-21) 


Here and in the rest of this part of Appendix, C denotes a generic positive constant independent 
of t and N. It follows from (IA-3al) . r n = v n and (|A-18I) that 


- 1 ElL\ " 1 I 

Pn{ rj Vn+ h 




®n+l) 


+ 4Ai 


Vn —1 V n 

r n -1 r n 


— ^ ('Pn+l ~ 'Pn) + C n , 


(A-22) 


where 


Vn = IPn 


r n - r n -1 


1/x n -l\ 2l (V n V n - l+2 V n . l 


r n -1 


r n - r n _i r n _i 


/ 

Vn 

1 

Vn Vfi — 1 

^71 Vfi — 1 

+ 

Vn 

\ 

( Pn 


Pn + T 

h 


-1- 



X n 

r n _i 

r n r n -1 


’Vn 

J 


(1A-21() follows from the summation of the product of r^_iV n and (1A-22|) . First, we analyze the 
second term on the left-hand side of (1A-22I) . Notice that 


TV—1 


TV-1 


n=l 


Y1 fa T n _]i) n — Y^ Y&n+l { r n V n +l r n-l^n) 

^ <3 +1 / ,„2 ^n+1 - 


7i=l 

N-l 


E ®» + i (+ AAi, | > 


71=1 


7V-1 


C EP 


71+1 


71=1 


/i h 

Vn+l V r 


n= 1 


4 ~ ‘JXn'Ur 


ex r . 


h 


+ {e + n l ) x n \v n 


and 




< Ce 


Vn+1 


/l 


A 


+ c 


Vn+1 V n 


h 


+ 


Xr, 
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It then yeilds from the Cauchy inequality that, for any 5 € (0,1), 


TV—l 


1 TV-1 ( 

X J (©n - ®n+l) r n-l^n > ~C5~ l £ N 


n=l 


+ 


TV-i ( , . 

EWPi 

n =1 l V 


TV-1 ( , 

2 Mpnr 

n=l l ^ 

Un. \ / .. 8. \ .o 4. \ u n+ i v n . 


~ V r . 


+ V„ 


+ ( 4A2 — -Ai \i> n + ( 4A2 + -Ai j x n -— 


TV-1 TV-1 ( , . . x 2 

C y (nT 1 + <5 _1 n~ 2 ) v 2 n - C{e + 25) y < x 2 ( n+1 J + 

n= 1 n=l l ' ' 


ft! 


(A-23) 


Similarly, 


TV— 1 / • 

1 | V n-1 _ ) r 2 


71—1 

^ r 

> - > <r. 


h \r n _ 1 r, ' ' n 

^71+1 ^71 


AT—1 . . AT—1 

E ^n—1 . V~^ ^Vi—1 T n T n — 1 .2 

r —i- h - + E T— 


E 

n=0 I 


h 


v n + hr n 


71—1 


f Vn+l 


71—1 


h Vn 


~ V n 


V /i 


TV-1 


+ y i^n - C (n 1 + e) u 2 } . (A-24) 


71—1 


So, if e is small (the smallness is independent of N) and N is large (which implies h is small), it 
follows from (IA-231) and (IA-241) that 


TV-1 


71 — 1 


— ®n+1 ) + ^Ai 

TV-1 


V n -1 V r 


+ ^E 


71—1 _ 


2 ( V n +1 V n \ .2 

I -7- I + v; 


V h 


r n -1 r r 

2 


TV—1 


r n-lVn > - 


c y -v 

^ n 


71 — 1 




(A-25) 


Here <r = min{2Ai/3, A2}. (Indeed, the hrst term on the second line of f|A-25[) is obtained by 
considering the following two cases: n = 1, • • ■ ,N — 2, and n = IV — 1, separately.) The first term 
on the right-hand side of (IA-251) can be bounded by 


TV-1 


[TV/2] 


TV-1 


[TV/2] 


TV—1 


TV-1 


I] ~fn < XX + jy X h 2 <C^p n h 2 + - J]! vl<y[Cp n + 


n= 1 n=l n=[TV/2]+l n=l n=[TV/2]+l 

This, together with (IA-251) . gives that for sufficiently large N, 


n=l 


AT/ Vn ' 


TV-1 

Eh 


71—1 


®n - ®n+l + 4A 


TV-1 


Vn-1 V r 


r n -1 r r 


N-l 


Vn-lVn > ~C y p n vl 


VE 


71—1 


2 f ^ 71+1 ^73 

V ft 


c£ K x 


x: 


n=l 
2 / ^n+1 V r 


TV-1 f /„. \ 2 


n=l f 


h 


+ Vr 


(A-26) 


In a similar but much easier way, we can deal with the other terms in (IA-221) and obtain (IA-211) . 
Step 2. In this step, we prove that 


TV-1 


hypnVn(t)+ / h 


71—1 


N / • / \ • / \ \ 2 N . 2 / \ 

VnW-Vn-lW^ , ^ 


E 

. 71 — 1 


h 


+ E 


n=l n 


<C£]\r(0) + C f [(^(s) + ^at( s )] 


(A-27) 
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Rewrite (IA-221) as 


/ \ 2 

Pn ( — ) V n A ~j^P (Qn Qn+ 1) = (P n+1 Pn) ~H Cnj 


(A-28) 


where 


^ _ ti n ti n —i | ^01—1 (^n V n ~i) 0 Vl 

Tn-rn-x r n _i (r n - r„_i ) 2 
I^nl E C (/5n|^n| I'^n/^-nl A |lVi/®n|) ■ 

Let 0 be a non-increasing function defined on [0,1] satisfying 
'0 = 1 on [0,1/4], 0 = 0 on [1/2,1] and |V’ / | < 32. 

Set 0 n = 0(x n ). Similarly, we consider the summation of the product of (I A-281) and tp n tin- To deal 
with the second term on the left-hand side of (j A-281) . one notices that 

N -1 JV-1 ^ 

' ~T (Qn Qn+1 ) 07101 — ' ~rQn (0 ntin 071—101—l) 


71=1 


71=1 

TV—1 


I i ti n ti n —i 0ti 0n—1 ■ 

Qn I 0n-7-1-J-^n-l | 5 


n=l 


h 


and 


TV — 1 . . . TV— 1 ^ / *2\ N—l . 

E Vn—l j V n ^71—1 V - ^ 1 ( I ^71^71 — 1 / V n \ V - ^ ^71 — 1 / V n ^71 — 1 

—^ —i — = 2^ h v n ~r r ~r~^ n+1 ~ J = ~ 2- — h — 

-j '71 — 1 IL i IL \ '71—1 '71/ -i '71—1 ,L 

71=1 71=1 x x 71=1 

N—l i / . • \ 2 TV—1 j i TV—1 I 

_ h t 71 / ^71-1 ~ ^71 \ t 71 — y71+l ^2 YU ^71 ~ ^71-1 ^2 


^ ^-1 V h 


71=1 

TV—1 


71=1 


r n -i hr n 


JV-1 . . . JV-1 , x . . v 2 JV-1 , , 

E ^n —1 / ^7i ^7i—l 7 Yn ( ^7i—l ^7i \ Yn Yn +1 -2 

j- 0”——s—J + L —— v » 

71=1 71=2 ' 71=1 


TV-1 

+ 2 0 

71=2 

which implies 


0n r n -r„_i . 2 
r i fir Un ’ 

' 71 — 1 71 


TV-1 . . . TV—1 . 2 TV—1 . / . . x 2 

V n -1 , Vn ~ Vn-1 ^ „/. r n ~ r n-1 < u I V n-1 ~ V r 

Yr\ 


n=1 


r n -i 


fi 


> 


71=2 

TV-1 

+ £ 

71=1 


"E 


h rl 1 r n _i 


71=2 


h 


0n ”071+1 . 2 

h^n 


Then, (IA-27D follows from (IA-211) . (j A- 6 [) . (IA-81) . (IA-91) and simple calculations. 
5fep 5. In this step, we prove 


JV-1 


£jv(t) Ah X] 1 


71=1 


^ti+i(^) 2u n (f) A u n _i(t) 


h 2 


A 


1 / u n 0 ) v n -x (t)' 


h \ X r 


%n—l 


<(721^(0) A C f [€;v(s) A 'S^-(s)] ds, t € [0, T\. 

Jo 


(A-29) 
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To this end, we rewrite (lA-3al) as 


P- 


Qn Gn+1 

~h 


— Pn 


V n + P 7 


exp{-7^ n+ i} - exp{-7t? n } 


+ 


-7 -7 

Pn +1 - Pn 


Here for n = 1, • • ■ ,N — 1, 
V n - r n _i 


<?n = In 


h 


+ 2 In 


1 - 


Ai—1 


r n+ i - ry 


h 

7 /„. \ 27' 


+ <?n 


- 1 =: 4 


|^n| is Cp n S |^n| T P 


7CI - 1 


^n—1 
1/ri £?n+l 


(A-30) 


(A-31) 


h 


+ 


r n +1 - r„ 


h 


- 1 


+ 


V rn 

—— -1 

%n. 


It is easy to derive from (IA-301) and (1A-27I) that 


TV—1 


, ^ | 9n(t) ~ Qn+\ (t) \ 

h ^Vn[ - h - ) 


< 


n= 1 

and then 

TV-1 


C<£n( 0)+C f [£jv(s) + Gi^r(s)] ds, 
Jo 


(A-32) 


h J2 Pn~ l (' gw ® <C£n(0) + cJ o [£ N (s) + £ 2 N (s)]ds. (A-33) 

Following the arguments in Section r3.4.1l (in particular, Lemmas 13.141 and 13.171) . we can use (1A-27I) . 
()A-32|) and fl A-.131) to get (|A-29|h 

Step 4■ Set T* = sup{t : (£j y(t) < 2(7(2bv(0)}, where C > 1 is the constant in (1A-34D . Obviously 
T* > 0 since C > 1. We may assume that T* < e/^2C<£ jv(0). (Otherwise, T* > e/\/4(7(5:(0) due 
to Cjv(O) < 2(£(0), and the lower bound is independent of N, so that the lemma is proved.) By the 
definition of T* , Gbv(i) < 2C'Gbv(0) for t E [0, T*]. Then, 

T* sup ^/£ n {s) < T*J2C£ n (0) < e. 
se[o,T*] v 

So, T* satisfies (1A-16I) . and it follows from estimate (IA-201) and definition of T* that 

2C€ n (0) = £ n (T*) < C£ n { 0 ) + C [ T [<£*(*) + <£^(s)] ds 

Jo 

< C£ n ( 0) + CT*2C<£ n (0 ) (1 + 2C , C JV (0)) . 

This implies 

T* > _L_ > _ 

“ 2C(1 + 2C£tv( 0)) - 2(7 (l + 4(72(0)) ’ 

The Lemma is proved by taking K = 2(7. 


□ 


Remark A-4 Let T* = min{T*, 1}. For the solution (r n ,v n ) obtained in Lem,ma \A-A it holds that 


TV-1 


£/v(t) + hJ2 P 2 n 1 


n =1 


^n+l(^) 2 V n (tf) T V n —i (t) 


h 2 


+ 


1 / U n (t) U n _l(f) 


h \ X r 


%n—1 


< ce(o) + ce 2 (o), t€[o,r*]. 

Indeed, (IA-341) follows from (IA-2911 and (1A-12D . 


(A-34) 
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Lemma A-5 Suppose that the assumptions in Lemma \A-2 I are satisfied. Then the following esti¬ 
mates are obtained for any n = 1, • • • ,N, and t,s € [0, X 1 *], 


1 < f n (t) < 3 1 < (r n - r n -i)(t ) < 3 


2 x n 
N 

*E 


n= 1 
AT 


“ 2 ’ 2 _ 


^n(i) - Tn-l(i) 


h 


h 

i 

+ h Yl 


- 2 ’ 


Vn(£) 




(u n - V n -l){t) 


N 'v n (t) - V n -l (t) 2 


*E 

72=1 

N-l 

72=1 

N-l 

fc E 


1 / r n (t) r n -i(t) 


h \ Xr 


%n—1 


72 — 1 
2 


/i 


< C7 (<£(0) + 1), 


< yci(o), 


+ 


^ 1 1 / w n (t) ■u n _i(t) N 2 


72=1 


h \ X r 


^72 — 1 


<c'x;« i (o), 

t=0 


1 /^_7-§ r n+l( i ) - r n (t) -i-\r n {t) - r n -i(t) 
Pn +1 i. Pn 


h 


h 


h 


< C((£(0) + 1), 


72=1 

N 


72=1 

AT 


/ ^n+l(*) ~ ^n(t) _ U n (2) ~ Vn-ljt) 

h \r n +i(t) - r n (t ) r n (f) - r n _i(t) 

N 

^\ v n(t) ~ •Un(s )! 2 + 

72=1 

r n [t) - r n -\{t) r n (s) - r n _i(s) 


N 'v n {t) v n {s) 2 


*E 

72=1 

N 

72=1 


3^72 ^72 

2 


i=0 

< C(£(0)|i — s|, 




h 


< C£(0 )|t-s| 2 , 


Vnit) - U n -l(i) fn(s) - ^ri-l(s) 


r n (t) - r n -i(t) r n (s ) - r n _i(s) 


< 


C^^(0)|t-s| +C£ 2 (0)\t-s\ 


(A-35) 

(A-36) 

(A-37) 

(A-38) 

(A-39) 

(A-40) 

(A-41) 

(A-42) 


2=1 


Here C is a constant independent of N. 


Proof. Clearly, (1A-35I) follows from (1A-12I) , (1A-14I) and (1A-18I) ; and (1A-36I) follows from (IA-121) . For 
(|A-37j) . it follows from (|A-35I) and (IA-341) that 


h T, 


N ' 1 / v n {t) v n -i (t) N 2 


72=1 

[JV/2] 


h V X r 


Xn —1 


N 


/[JV/2] 

* E + E I 

\n=l n=[N/2]+lJ 


) 


1 / V n (t) V n -l(t) 


h V X r , 


%7l — l 


<Ch £ 


- 27 -i 


72=1 


1 / w n (i) 


h V Xr, 


^72 — 1 


AT 


+ Cft XI 

n=[JV/2]+l 


Vnit) - V n -i (t) 


h 


+ 


Vnit) 


<C£(0) +C£ 2 (0). 

Similarly, the estimates for r n in (IA-3711 follows from (I A-35[) and (IA-121) . (I A-38[i follows from simple 
calculations and (IA-121) . For ()A-39I) . we note that 

1 V n+ i(t) - Unit) Vnit) - Vn-lit) 


r n +i it) - r n it) r n (t) - r n _i (t) 

< 


g n+1 - 2^) - ( g n - 2^ 

r n J \ r n -1 


Gn+1 Gn 

♦! 

( X n 

V n 

V n -1 

+ 

Vn-l 

%n %n—l 

r n 

r n -1 

h 

K r n 


^72—1 

%n—l 

r n r n ~i 


3^72 — 1 


Then (IA-391) follows from (IA-321) . (|A-35|) and (1A-37D . For (I A-40 p . notice that 


N 


*E 

72=1 


Vnit) V n (s) 


x r . 


N 


h T, 

72=1 


Vnir) 


dr 


Xr 


JV rt 

s*E/ 

n=l Js 


Vnir) 


dr\t — s|. 
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Then the estimate for v n /x n in (IA-401) follows from (1A-27I) . With this, the estimate for v n in (IA-401) 
holds obviously. Similarly, (|A-41|) follows from (1A-35D : (1A-42I) from (lA-27jl and (1A-35D . 

□ 


For h = 1 /N, we dehne the functions (r h , v h ){x, t) as follows: 

h( + \ , + \ , r n{t) ~ r n -i(t) , . 

r n (x,t ) = r n _i(t) H--- (x - x„_i), 

v h (x,t) = v n -!(t) + ln ^ ^ n ~ l ^\ x - X n -l), 

for x n -\ < x < x n , 1 < n < N and 0 < t < T*. Then we have 


_ K, .s fc/ _ r n (t) - r n -i(t) u, _ v n {t) - V n -! (t) 


r?(x,t) = v h (x,t), r%(x,t) = 


h 


and v^(x,t ) = 


h 


Proposition A-6 If the assumptions in Lemma \A -HI are satisfied, then there exist subsequences of 
{r h }, {v h }, {r h /x}, {v h /x}, {/0 7 ~ 1//2 r^} and {v^/r^:}, still labeled by {r h }, {v h }, {r h /x}, {v h /x}, 
{p7-!/2 r ft.} anc i {v!f/r!f} for convenience, such that {r h }, {v h }, {r h /x}, {v h /x}, { p 7 _1 ^ 2 r^} and 
i v x/ r x } converge boundedly and almost everywhere in [0,1] x [0,T*] for h —> 0. 


Proof. We use similar arguments as in p32l 25] . First, we consider {r h } and It follows 

from (|A-35jl and ()A-36I) that the functions of the families {r h } and as functions of x have 
uniformly bounded total variations with respect to h for each fixed time t € [0,T*]. Let t = tk 
[k = 1, 2, ■ ■ •) be a countable set which is everywhere dense in [0,T*]. By Helly’s theorem and 
a diagonal process, from the family of functions {v h }, we can select a subsequence, still labeled 
as for convenience, converging boundedly and almost everywhere in x € I on the dense set 
{tk',k = 1,2, • • • } in [0,T*] as h —>• 0. Consequently, by the Lebesgue’s theorem, the subsequence 
{u^} converges in L 2 -norm on {tk', k = 1, 2, • • • }. Next, by (| A-40[) . the continuity in time for the 
L 2 -norm of v h , it is standard to show that the subsequence converges in L 2 (I) uniformly 

in t £ [0, T*], as h -> 0. So, we can select further a subsequence, again still labeled as {u^ 1 } for 
convenience, converges almost everywhere in (x,t) € I x [0, T*]. Denoting the limiting function of 
{u^} by v. Since rjf(x,t) = v h (x,t), we know that the corresponding subsequence of {r h } (still 
labeled as {r h } for convenience) converges almost everywhere to r(x,t ) := ro(x) + f ( * v(x, s)ds in 
/ x [0, T*]. 

Similarly, one can derive from estimates (|A-35|) . (IA-371) and ()A-40H that there exist certain 
subsequences of {r h /x} and {v h /x}, still labeled as {r h /x} and {v h /x} for convenience, converge 
boundedly and almost everywhere in (x,t) £ I x [0, T*} to the functions r(x,t)/x and v(x,t)/x, 
respectively. The estimates flA-351) . (IA-38D . (IA-391) . (IA-41D and (IA-42D guarantee that one can 
select certain subsequences of {p 7-1 ^ 2 ^} an d { v x/ r x\ (still labeled as {p 7 ” 1 ^ 2 ^} anc ^ { v x/ r x} f° r 
convenience) such that they converge boundedly and almost everywhere in (x, t) £ I x [0, Tfi to 
the functions fP~ 1 ^ 2 r x (x,t) and v x (x,t)/r x (x,t), respectively. 

□ 

Due to the above convergence and uniform estimates on the approximating sequence, one may 
first verify that the limiting function (r, v) is a weak solution to problem (12.1.1111 in the sense of 



4vr y z p(y)dy 


(x, t)dxdt 


— + 2-) if x + 2p(-Y -vA (x,t)dxdt 
r x r J Vr/xJ 
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for any test function ^ £ C 1 ([0,T*] x I) satisfying £ C c (I). Then one can use the standard 

regularity argument (cf. [26]) which is ensured by the above uniform estimates to show it is also a 
strong solution as in Definition 12.11 satisfying 

l£(t) £ C([0,T*]) and f —— 03^ (x, t)dx < C£(0), t£[0, T*]. 

Ji P{x) 

This in particular implies that ©(•,£) £ H l (I ) for t £ [0,T*], so that one can define the trace of 
at x = 1 . 

The uniqueness of the strong solution follows from the argument in |29] (Section 11) of the 
weighted energy estimates. 


Part II. Linearized analysis 


We use the simple case that A 2 = (2/3)Ai = 1/3 to illustrate the main ideas of the linear analysis. 
As in na. we set w = r/x — 1. Then 

r = x + xw, r x = 1 + w + xw x , v = xw t , v x = w t + xw tx , v t = xw tt . (A-43) 

The linearized problem for (12.1. lip around w = 0 (the equilibrium of (12.1.111) ) reads 

xpwu - (3y - 4 ){p 1 ) x w - 7 a ; -3 ( p y x 4l w x ) x = (xw xt + 3 w t ) x , (A-44a) 

53^ := xw xt + w t = 0 at x = 1. (A-44b) 


The condition v(0,t) = 0 has be incorporated in the transformation from r to w since v = xwt- 
(Indeed, (IA-441) follows from the facts that the principal parts of ( r/x)~ 2j , r)) 7 , (r/x)" 4 , v x /r x and 
v/r are 1 — 2^w, 1 — 7(10 + xw x ), 1 — 4 w, xw x t + wt and wt , respectively. The derivation of the 
left-hand side of (|A-44aP can be also found in HZ]-) Naturally, we may rewrite (|A-44al) as 

x A pwu + (3y - 4 )4>x 4 pw - 7 (x 4 ^ 7 ^)^ = x 3 (%$l) x + 2x 3 wtx, (A-45) 

where (f> is dehned in ( 12 . 1 . 2 p . which is bounded from above and below by positive constants, so that 
it is easy to see that 7 > 4/3 is crucial to the linearized stability. 

Lower-order estimates. The basic multipliers for f)A-45[) are wt and w which yield the bound¬ 
edness of 


(r - x, xr x - x)(•,t )|| 2 + (1 +t) \\(v,xv x ) (-,t )|| 2 + (1 + t) xp 1/2 v t {-,t) 


(A-46) 


Indeed, we have the following basic estimates for the solution of (|A-44D if 7 > 4/3. 

J x 4 [pw 2 + (37 - 4 )c/)pw 2 + 7 p^w 2 ) dx + J x 2 [(w t + xw tx ) 2 + 2 w 2 \ dx = 0 , 

J [x 2 (w + xw x ) 2 + 2x 2 w 2 + x A pwwt\ dx + J [(37 — 4)0x 4 pu; 2 + r yx A p 1 w‘^\ dx 


which implies 
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Since the coefficients for the equation and boundary conditions in (IA-44D are independent of t, 
differentiating (1A-44I) with respect to t, one obtains the same estimates for the corresponding 
t-derivatives of each quantity appearing in (|A-47H . in particular, 

(1 + t) J x 4 (pwu + P w t + f w tx) ( x i t)d x + J (1 + s) J x 2 [(io ss + xw ssx ) 2 + WsJ dxds 
< C j [x 4 w 2 + x 2 w 2 + x 4 pw 2 + x 4 pw 2 t + x 4 p ) w 2 x \ (x, 0 )dx. (A-48) 

Integrating the identity 

^ (1 + t) J x 2 [( w t + xw tx ) 2 + w 2 ] (. x , t)dx = J x 2 [(w t + xw tx ) 2 + w 2 ] (x, t)dx 

+2(1 +t)Jx 2 \(wt + »„)(»» + xwttx) + «-.»«-] (+, *)<*+ 
with respect to t, and use (lA-47ji and (1A-48I) to get 
(1 + t) J [x 2 {w t + xw tx ) 2 + x 2 w 2 ] (x, t)dx 
<C J [ x 4 w 2 + x 2 w 2 + x 2 w 2 + x 4 pw 2 t + x 4 w 2 x ] (x, 0 )dx. 


(A-49) 


Therefore, the boundedness of (1 A-55[) is a consequence of (|A-47H - (|A-49[) and (IA-431) . 
Higher-order estimates. We may rewrite (|A-44al) as 


G x t + "fP^Gx = xpw u + 'y4>x 2 pw x + (3y — 4 )<f>xpw, where G = xw x + 3 w. 


(A-50) 


Based on (j A-47h . (IA-481) and ()3.2.9I) . we may apply a multiplier p a G x with a > 2y — 2 to (1A-50I) to 
get 


(A-51) 


J p a G 2 (x,t)dx + J J p^ +a G 2 x dxds 

<C J [p a G 2 + x 4 w 2 + x 2 w 2 + x 4 pw 2 + x 4 pw 2 t + x 4 fPw 2 x ] (x, 0 )dx. 

Due to equation (IA-501) . we have 

J x 4 pw 2 t (x,0)dx — Gj x 2 p~ l (\G x t\ 2 + Ip'Gx] 2 + \x 2 pw x \ " + \xpw\ 2 ^j (x,0)dx, 

which implies that the largest a in (IA-51D could be 2 y — 1 . So, we choose a = 2 y — 1 . Moreover, it 
follows from (1A-51D and equation ()A-50|) that 

J J p^~ l G 2 xt dxds — G J [p 2 ' l ~ 1 G 2 + x 4 w 2 + x 2 w 2 + x 4 pw 2 + x 4 pw 2 t + x 4 p' y w 2 x ] (.x , 0 )dx. 
This may suggest 

f p 2 ' y ~ 1 (x 2 w 2 x + w 2 )(x,t)dx + f f p 3 ' y ~ 1 (x 2 w sxx + w 2 sx )dxds 

J Jo J (A-52) 

<C J [ p 2l ~ l G 2 + x 4 w 2 + x 2 w 2 + x 4 pw 2 + x 4 pw 2 t + x 4 p' y w 2 x \ (x, 0 )dx. 
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(Indeed, it is one of main ideas to justify this for the corresponding nonlinear equation). The 
higher-order estimate (|A-52jl can improve the regularity near the origin as follows. It follows from 
(1X321) . (1X471) and (EHl that 


f (w 2 + w 2 xs )dxds <cj [p 27 1 G 2 x + x 2 w 2 + x 4 ( w 2 + pw 2 t + pwf t + f ! w 2 tx )] (x, 0 )dx. (A-53) 

Let ^ be a non-increasing cut-off function defined on [0,1] satisfying if = 1 on [0,1/4], ip = 0 on 
[1/2,1], and \ip'\ < 32 . With (|A-53[) and (1A-48|) . we can integrate the product of (IA-44al) t and 
ipxwtt with respect to x and t to get 

/*!/ 4 r 

J x 2 w 2 t (x, t)dx < C J [p 27 ^ 1 ^? 2 + x 4 w 2 + x 2 w 2 + x 4 pw 2 + x 2 pw 2 t + x 4 p y w 2 x \ (, x , 0 )dx, 
which, together with (1A-48I) . gives 

J x 2 pw 2 t (x, t)dx < C j [ p 2l ~ l G 2 + x 4 w^ + x 2 w 2 + x 4 pw 2 + x 2 pw 2 t + x 4 p' y w 2 x \ (x, 0 )dx. (A-54) 
As a consequence of (I A-43[) . (|A-52j) and (|A-54[) . we can obtain the boundedness of 


P 1 / 2 vt(-,t) 


+ 


57-V 2 


(r xx , {r/x) x )(-,t) 


(A-55) 


Conclusion: A natural higher-order functional for the study of the linear problem (IA-44D is: 


E(t) = || (r-x,xr x - a;)(-,t )|| 2 + \\{v,xv x ) (-,t )|| 2 + p 1/2 v t (-,t ) 


+ 


P 7 1/2 {r xx , (r/x) x ){-,t) 
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